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AN OPTIMIZED PRINTER PLOTTING SYSTEM CONSISTING OP 


COMPLEMENTARY 7090 (FORTRAN) AND 
1401 (SPS) SUBROUTINES 
PART II - SYSTEMS PROGRAMMER MANUAL 
by Lois T. Dellner and Betty Jo Moore 
Lewis Research Center 


SUMMARY 

Complementary subroutines, easy for the FORTRAN programmer to use, provide 
printed plots as part of the normal off-line output. These subroutines are simi- 
plified so that no choices must be made and generalized so that choices may be 
made, by the programmer, of the plotting characters, the scales, the appearance 
of the grid, the number of curves, and other options. The 7090 routines write 
information-packed records that axe decoded and used by the 1401 subroutine to 
produce the plot at full printer speed (6 sec/page) . Part I (Technical Note 
D-2174) of this presentation is an instruction manual for the users of the sys- 
tem. Part II is intended to simplify the work of the systems programmer who im- 
plements, and probably must modify, the system, and then must help debug users’ 
programs . 


INTRODUCTION 

v The 0PP system (Optimized Printer Plotting), a set of complementary sub- 
routines for the IBM 7090 and 1401, is a simple-to-use system that requires mini- 
mum storage, reduces 7090 execution time (by a factor of more than 10 from that 
required by presently available machine plotting methods), and produces printed 
plots at full printer speed (6 sec/page) as part of normal off-line output. This 
optimization is achieved by combining three techniques: (l) performing most of 

the requisite arithmetic on the 7090, (2) packing the results on the output tape, 
(3) having the 1401 generate the grid and the grid-labels, as well as position 
the points in a subroutine of the standard tape-to-pr inter program. 

The system permits, but does not require, the programmer to choose the 
plotting characters, the scales, the grid-line spacing, etc. It handles single 
or multiple curves, prints true scales, and permits printing of titles above and 
below the plot. 

The chief advantages of 0PP system plotting over manual plotting are the 
complete elimination of elapsed time and the reduction of the cost per plot by a 
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factor of at least 50. Even when the limitations of the printer - 10 characters 
per inch horizontally (to a minimum of 100) and 6 characters per inch verti- 
cally - make a printed plot inadequate for a particular application, the time 
required to produce a manual plot from the same data is reduced by at least 
5 minutes, because the printed plot supplies guidance in the choice of scales 
and the minimum number of points required to define the curve. 

More than 150 programmers are now using this system effectively and more 
than 500 plots are printed weekly. The complementary subroutine PLOT for the 
1401 was written for a 4K machine. 


GENERAL DESCRIPTION 

The 0PP system offers the programmer plotted' output with a minimum of ef- 
fort. After writing a title for his plot on the output tape, he writes CALL 
PL0TXY or CALL PL0TMT (for multiple curves). The arguments, or call list, in- 
clude the names of the arrays to be plotted and specify the number of points per 
curve and the number of curves. The programmer then writes a legend to be 
printed at the bottom of the plot. 

The plot or plots are printed as part of the regular output listing with no 
delay. No changes in his Card-handling procedures nor special instructions for 
the operators are required. 

If he is using PL0TXY, the values of the variable to be plotted in the 
x -direction must be in sequence. If they are not, the subroutine S0RTXY is sup- 
plied to be used before calling PL0TXX. For either PL0TXY or PL0TMT, if the 
size of the elements in (or the total range of) any array is not known to be 
within certain limits, the programmer calls the subroutine SCALE for each array 
before calling the plotting subroutine. SCALE will transform the array to suit 
PL0TXT and FL0TMY only if it is necessary. 

In addition to the minimum-effort use just described, the programmer may 
choose to use one or more of several options that permit him to control, for 
example, the appearance of the grid (by specifying the frequency of the grid- 
lines In either direction) , the scale for either variable (by specifying the 
scale-factor and a starting-value), the plotting character, etc. 

This presentation of the 0PP system is in two sections. Part I (ref. l) 
is a manual for the FORTRAN programmer who wants printed plots as part of his 
normal output. The general objective of Part II is to simplify the work of the 
systems programmer who Implements, and probably must modify, the 0PP system, and 
then must help debug users' programs. To such a programmer, the system consists 
of a set of FORTRAN II subroutines (PL0TXY, PL0TMT, PISTUG, SCALE, and S0RTXY) 
and one 1401 (SPS) subroutine (PLOT) incorporated into a standard tape-to-printer 
program. PLOTXY and PLOTMY write information-packed records on the output tape, 
and PLOT decodes these records and generates the plot, one line at a time, pre- 
paring the next line as the current line is being printed, maintaining full 
printer speed the while. 


2 



A brief discussion and a complete listing are included for all routines* 

For IhOTXY, PLOTMY, and PLOT, detailed specifications of input and output and 
comprehensive block diagrams are also provided* For PLOTXY, PLOTMY, and PISTUG, 
dictionaries of the FORTRAN variables are supplied* These are intended as sys- 
tems debugging aids* Each is designed for use primarily as a cross-reference 
with the corresponding FORTRAN listing and block diagram. Within each dictio- 
nary, each variable name is identified in terms of its contribution to the pro- 
gram. Throughout the text, a subroutine name is printed at the bottom of each 
page for assistance in rapid cross-referencing. 

Appendix A consists of a set of sample plots each of which is accompanied by 
a representation of the associated tape records written by the 7090 and used by 
the 1401 to produce it* The sources of some difficulties are discussed in ap- 
pendix B, and corrective actions suggested where feasible. 
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PLOTXY - (FORTRAN II) 


To plot one curve, the programmer will usually use PLOTXY. His call state- 
ment lists the names of the arrays to he plotted, the number KODE, (an indicator 
of which options are being used (see I.-D.), and finally, an array whose first 
element is the number of points to be plotted and whose other elements are any 
additional data required by the choice of option, 

PLOTXY may call an auxiliary subroutine, PISTUG, (p. 44) to compute scaling 
parameters, then it writes one INFO record specifying the appearance of the grid 
to the 1401 subroutine PLOT. This is followed by one or more packed DATA records 
that specify the location of each point to be plotted in terms of KX, the number 
of lines down from the first x grid-line, and KY, the number of print positions 
to the right of the first y grid-line. 

I. Input to this routine is supplied via the four arguments of the CALL, as 
follows : 


CALL PLOTXY (XDOWN, YACROS, KODE, P) 

A. XDOVJN and YACROS are floating-point variable array names. In general, 
the values in both arrays must lie between -10® and 10®. (See appen- 
dix B, VI., for precise limitations,) The subroutine SCALE (p. 54) is 
supplied to handle arrays whose element size Is not known- In this case, 
SCALE should be called by the programmer before calling PLOTXY. 

The elements of the XDOOT array must be in either decreasing or 
increasing sequence. The subroutine SORTXY (p. 52) is supplied to 
rearrange the arrays when the XDOOT array Is not in order. It should 
be called by the programmer before calling PLOTXY. 

B. KODE is a fixed-point variable or constant whose value is restricted 
by 0 < KODE < 127. Seven different options are offered to the pro- 
grammer. Each has a number associated with It. The sum of the numbers 
corresponding to the options used is the KODE. 

C. P is a floating-point variable array name. P(l), the first value in 
the P array, is the number of points to be plotted. The other values to 
be supplied in the P array are determined by the particular options being 
used. 
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D. Options. The sample plots in appendix A illustrate the effects 
of using the various options. 


NUMBER 


2 


4 


16 


OPTION 

Choic e of plotting character . - The plotting character will 
be an asterisk unless the programmer uses this option. When 
he does, he must specify the plotting character by placing 
its BCD code in the leftmost six bits of P(2) . He may use 
any character in the set of all FORTRAN characters except 
the minus sign ,r -". 

Choice of x grid-line frequency . - An x grid-line will be 
printed ^ very 10 line-spaces below the first one unless the 
programmer uses this option. When he does, the desired 
frequency (M) of x grid-lines must be supplied to PLOTXY in 
P(3). If M » 0, only the first and final x grid-lines will 
be printed. 

Choice of y grid-line frequency . -Ay grid-line will be 
printed every 10 print positions to the right of the first 
(leftmost), -unless this option is used. If it is, the de- 
sired frequency (N) of y grid-lines must be supplied in 
P(4). If N = 0, only the first y grid-line will be printed. 

Choice of x scale . - PLOTXY will supply the three scaling 
parameters that specify the x scale star ting -value, the 
scale-factor (value of one line space) , and the x grid- 
labels, unless this option is used. When it is used, the 
programmer must specify his own scale by supplying the three 
parameters KSX, FX, and DX in P(6), P(7), and P(8), respec- 
tively. These parameters are defined by: 

FX Is the desired starting -value multiplied by 10*^ . 

DX is the desired scale-factor for one line-space 
multiplied by 10^. 

U must be chosen so that FX and DX are integers and 
0 < U < 6. 

KSX is equal to 6 - U. 
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NUMBER 


OPTION 


32 Choice of y scale . - PLOTXY will supply the three scaling 

parameters that specify the y-scale starting -value, the 
scale-factor (value of one print position), and the y grid- 
labels, unless the programmer uses this option* If he /does, 
he must specify his own scale by supplying the three param- 
eters, KSY, FY, and DY in P(9), P(lO), and P(ll), respec- 
tively- These parameters are defined by: 

FY is the desired starting -value multiplied by 10^. 

DY is the desired scale-factor for one print position 

multiplied by 10^. 

V must be chosen so FY and DY are integers and 

0 < V < 6. 

KSY is equal to 6 - V. 

64 Choice of printing coordinates at left of plot * - The 29 

print positions to the left of the first y grid-line are 
usually blank except for x grid-labels every tenth line in 
positions 20 to 28- If option 64 Is used, the coordinates 
of each point are printed immediately to the left-*- of the 
grid-labels on the same line as the point* 

8 Choice of printing message at left of plot * - The 29 print 

positions to the left of the first y grid-line are usually 
blank except for x grid-labels printed in positions 20 to 28. 
When option 8 is used, the programmer must supply a 6- in 
P(5) and the information to be printed in P(l2), P(l3), 

.... The 6 BCD characters in P(l2) are printed immedi- 
ately to the left- 1 - of the grid-labels on the same line on 
which the first point is printed, the six characters in 
P(l3) are printed with the second point, etc- Specifically, 
the 6 BCD characters in P(l+ll) are printed on the same line 
on which the I^ 1 value of the XDOWN array is plotted. 


^If options 8 and 64 are used together, the BCD characters will be printed 
to the left of the coordinates. 
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II. Output from this routine consists of one INFO record and one or more 
DATA records. Sample tape records appear in appendix A. 


A. INFO record. 

The contents of the INFO record are the information needed by the 
1401 subroutine PLOT to determine the construction of the grid, the con- 
struction of the grid-labels, the plotting character, and the format of 
the succeeding DATA records. Specifically, the INFO record consists of 
the two characters PI, followed by 10 fields, whose contents are: 

(Field width is shown in parentheses.) 

a. (l) This is the plotting character for all points in the DATA 
records . 

b. (l) This digit is the decimal -point-shift code, KSX, for OFX and 

ODX. (items d. and f. in this list.) 

c. (l) This digit is the decimal -point -shift code, KSY, for OFY and 

ODY. (Items e. and g. in this list.) 

d. (9) The number (OFX) in this field is in the form iO.xxxxxx and 
represents the x-scale starting -value. 

e. (9) The number (OFY) in this field is in the form IO.xxxxxx and 
represents the y-scale starting -value. 

f. (9) The number (ODX) in this field is in the form IO.xxxxxx and 
represents the x scale-factor (value of one line-space). 

g. (9) The number (ODY) in this field is in the form IO.xxxxxx and 
represents the y scale-factor (value of one print position) . 

h. (3) This integer is the frequency, M, of the x grid-lines. 

i. (3) This integer is the frequency, N, of the y grid-lines. 

j. (3) This integer specifies the length, NF, of the first field in 
each block of information in the following DATA record. NF may be 
zero. 
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B. DATA records. 


The contents of a DATA record (after the first two characters PD) 
are successive blocks of the information needed by the 1401 subroutine 
PLOT to position one plotting character and print the left-side field on 
the same line. Within a record, only complete blocks are permitted; one 
record may hold as many as 21 blocks. Within a block, the location in 
which the plotting character for one point is to be printed is specified 
by two positive integers, KX (the number of line-spaces below the first 
x grid-line) and KY (the number of print positions to the right of the 
first y grid-line) . Specifically, the format of one block in the DATA 
record is: 

a. The first HP (see j. preceding) characters are those to be printed 
in the left-side field (to the left of the first y grid-line) unless 
HP is zero; in this case, the characters to be printed are blanks and 
will be generated by the 1401 subroutine. 

b. The next three characters are KX. 

c. The next (and last) three characters are KY. 

III. After "initializing" switches, the routine inspects the input variable 
KODE bit by binary bit, using Boolean statements and octal masks. It 
branches on each bit to set necessary output switches, and either supplies 
standard values or picks up programmer -provided specifications from the P 
array. If necessary. It calls an auxiliary routine (PISTUG) for scaling 
either or both coordinate arrays. The INFO record (II.-A. ) is then written. 

The routine then uses the scaling parameters to calculate KX (the num- 
ber of line-spaces down from the first x grid-line) and KY (the number of 
print positions to the right of the first y grid-line) for each point, which 
together determine where the plotting character should be placed. As soon 
as enough points have been processed to fill a DATA record, it is written 
out, and this continues until all the points have been processed. Control 
is then returned to the calling program. 
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SUBROUTINE PLOTXY ( X , Y , BK , P ) 


to 


c 

c 


c 

c 


c 


c 


c 


c 

c 

c 

B 


c 

B 


C 

B 


C 

B 


C 
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SUBROUTINE PLOTXY(X,Y,BK,P) 


THIRD ARGUMENT IS A BOOLEAN VARIABLE 


DIMENSION X{ 100) f Y( 100 ) , P (1 11 ) , KX( 21 ) ,KY ( 21 ) , XOW< 10) , YOW(IO), 
1 PQW(IO) 


EQUIVALENCE ( BKK ,KK ) 


IT IS ALSO USED AS AN INTEGER 


100 N*P( 1 ) +SIGNFI *5* P( 1 ) ) 

101 BKK=BK 

102 KSW=0 

103 KHDG=2 

GO TO 140 IF KODE IS NOT ZERO 

104 IF(KK) 140,106,140 

106 A 1 = 1H» 

107 M= 10 

108 NY = 10 

109 NF 1=0 

SCALE X AND TO 190 OR ERROR EXIT 

120 K=3 

121 CALL PISTUG(K,N,X,DX,FX,KSX,KHARX) 

122 KDX=DX 

124 IF(KDX) 190,900,190 

SCALE Y AND TO 220 OR ERROR EXIT 

126 K=2 

128 CALL P1STUG(K,N, Y,UY,FY,KSY,KHARY) 

130 KDY=DY 

132 IFUDY) 220,900,220 


OPTION 1 - USE • OR P ( 2 ) 

140 I F ( BKK*1000000 ) 148,148,144 
144 A 1-P ( 2 ) 

146 GO .TO 150 
148 A1=1H* 

OPTION 2 - USE 10 OR P ( 3 ) 

150 I F ( BKK*2000000) 158,158,154 
154 M=P(3)+SIGNFl.5,P{3)) 

156 GO TO 160 
158 M=1 0 

OPTION 4 - USE 10 OR P(4) 

160 IF ( BKK*40 00000 i 168,168,164 
164 NY = P 1 4 ) +S IGNF t . 5, P ( 4) ) 

166 GO TO 170 
168 NY * 10 

OPTION 8 - ALTER KSW AND NF1 

170 I F( BKK*1QQOOQOO ) 178,178,174 
174 NF1=P(5) + SIGNF( • 5, P< 5 ) ) 

176 KS W= 1 

177 GO TO 200 

178 NF 1=0 

OPTION 64 - ALTER KSW AND NFl 

200 IF ( BKK* 100000000 ) 180,180,204 
204 NF 1=NF 1+13 
206 KSW=KS W+2 
208 KHDG=1 
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SUBROUTINE PLOTXY ( X, Y, BK , P ) 


C OPTION 16 - GO TO 120 OR USE P ( 6-8 ) 

B 180 IF(BKK*20000000) 120,120,184 

184 KSX=P16)+$IGNF( • 5, P ( 6 ) ) 

185 F X = P ( 7 ) 

186 DX = P { 8 ) 

187 RFX=FX*( 10.**(KSX-6) ) 

RQX=DX* ( 10 • ** ( KSX-6 ) ) 

PFX=RFX+100.*RDX 

XMAX=MAX1F{ AB$F( RFX) ,A8SF(PFX) ) 

KHARX=XlNTF(LOGF( XMAXJ/2. 302585+40. ) -40 

C OPTION 32 - GO TO 126 OR USE P(9-ll) 

B 190 IFl BKK+40000000) 126,126,194 
194 KSY=P(9)+$IGNF( .5,P(9) ) 

196 F Y = P ( 10 ) 

198 DY=P ( 11) 

199 RFY=FY*( 10.**tKSY-6) ) 

RDY=DY*( 10.**(KSY-6) ) 

PFY=RFY+10Q.*RDY 

YMAX=MAX1F( ABSF(RFY) , ABSF ( PFY ) ) 

KHARY=XINTF(LOGF( YMAX)/ 2. 302585+40. ) -40 

C GO TO 250 OR PRINT SPECIAL HEADING 

220 GO TO ( 230,250) ,KHDG 
230 KPWRX=3-KHARX 
232 KP WRY=i-KHARY 
234 F10X=IQ.**KPWRX 
236 F10Y=10.**KPWRY 

238 WRITE OUTPUT TAPE 6, 500 , KPWRX , KPWR Y 

250 ODX = DX* l.E— 6+SIGNF( .0000005, DX ) 

OFX = FX* 1 • E— 6 + S IGNF( .0000005, FX ) 

ODY = UY*1 .E-6+SIGNF(. 0000005, DY ) 

OF Y = FY* 1 • E-6+S IGNF ( .0000005, FY) 

C WRITE INFO RECORD 

251 WRITE UUTPUT TAPE 6, 502, A1 , KSX, KSY , OFX , OFY , ODX, ODY,M, NY,NFi 
C 
C 
C 

300 L= 1 30/ ( 6+NF 1 ) 

J=1 
C 

310 DO 498 1=1, N 
C 

320 C6 = XU)*t 10.**(6-KSXi i/DX 
322 C7=INTF(C6+SIGNF( • 5 , C6) ) 

324 KX(J)=(DX*C7-FX)/DX+.5 
330 C6 = Y ( I ) * ( 10.**(6— K$Y) )/DY 
332 C7= INTF ( C6+S IGNF ( . 5 * C6 ) ) 

334 KYI J>=lQY*C7-FY)/DY+.5 
C 

338 KBR=KS W/2+ 1 
340 GO TO (346,342) ,KBR 
342 XOWl J)=F10X*X( I) 

344 YOW ( J ) =F10Y*Y ( I ) 

C 

346 KBR2=KSW+1 

347 GO TO(350,348,350,348),KBR2 

348 POW( J)=P( I+li) 


INITIALIZE FOR OUTPUT LOOPING 
FOR ALL POINTS 

COMPUTE COUNT DOWN AND ACROSS 


GO TO 346 OR PICKUP COORDINATES 
GO TO 350 OR PICKUP P ARRAY MESSAGE 
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SUBROUTINE PLOTXY { X , Y ,BK , P ) 


C 

C 

C 


c 

c 


GO TO 360 UNLESS READY TO WRITE RECORD 

350 I F ( J-L ) 360,352,360 

352 J=0 

354 GO TO 400 

GO TO 498 IF ALL POINTS NOT DONE 

360 I F ( t-N ) 498,362,498 
362 L=J 

WRITE A OATA RECORO 


400 KBRM^KbW+i 

402 GO TO 1410,420,430, 440) , K8RM 

410 WRITE OUTPUT TAPE 6 , 510 , ( KX ( J ) t KY ( J > , J=1 , L ) 

412 GO TO 498 

420 WRITE OUTPUT TAPE 6,520, ( POW( J ),KX(J),KY(J),J=1,L) 

422 GO TO 498 

430 WRITE UUTPUT TAPE 6, 530 , ( XOW (J > , YO W( J ) , KX (J ) , KY ( J ) , J=1 , L ) 

432 GO TU 498 

440 WRITE UUTPUT TAPE 6, 540 , ( POW t J ) , XO W(J ) , YOW ( J ) , KX ( J ) , KY l J ) , J=1 , U 
498 J=J+l 

END OF DO THAT STARTED AT 310 

900 RETURN 


500 FORMAT (2HPT/2HPT/2HPT,8X,3HX*E» 12, 5H Y*E,I2) 
502 FORMAT 1 2HP I ,A1 , 2 1 1 , 4F9.6, 313 1 
510 FORMAT (2HP0,42I3> 

520 FORMAT t 2HP0, 101 A6 , 13, 13) > 

530 FORMAT t 2HP0, 6{ F6.0,F7.0, 13,13)) 

540 FORMAT ( 2HP0 , 5 ( A6 , F6. 0 , F7. 0 , 13,13)) 

END 11,0,0,0, 0,0, 0,0, 0,1, 0,0, 0,0,0) 


H 

H 
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SUBROUTINE PLOTXY ( X, Y* Bk , P ) 


STORAGE NOT USED BY PROGRAM 


DEC OCT DEC OCT 

875 01553 32561 77461 


STORAGE LOCATIONS FOR VARIABLES APPEARING IN DIMENSION AND EQUIVALENCE STATEMENTS 



DEC 

OCT 


DEC 

OCT 


OEC 

OCT 


OEC 

OCT 


DEC 

OCT 

BKK 

874 

01552 

KK 

874 

01552 

KX 

873 

01551 

KY 

852 

01524 

POW 

811 

01453 

XOW 

831 

01477 

YOW 

821 

01465 











STORAGE LOCATIONS FOR 

VARIABLES NOT APPEARING 

IN COMMON* DIMENSION* OR EQUIVALENCE 

STATEMENT 




DEC 

OCT 


DEC 

OCT 


DEC 

OCT 


DEC 

OCT 


DEC 

OCT 

Ai 

801 

01441 

C6 

800 

01440 

C7 

799 

01437 

DX 

798 

01436 

DY 

797 

01435 

F10X 

796 

01434 

F10Y 

795 

01433 

FX 

794 

01432 

FY 

793 

01431 

I 

792 

01430 

J 

791 

01427 

KBR2 

790 

01426 

K8RM 

789 

01425 

KBR 

788 

01424 

KDX 

787 

01423 

KDY 

786 

01422 

KHARX 

785 

01421 

KHARY 

784 

01420 

KHDG 

783 

01417 

KPWRX 

782 

01416 

KPWRY 

781 

01415 

K 

780 

01414 

KSW 

779 

01413 

KSX 

778 

01412 

KSY 

777 

01411 

L 

776 

01410 

M 

775 

01407 

NF1 

774 

01406 

N 

773 

01405 

NY 

772 

01404 

ODX 

771 

01403 

ODY 

770 

01402 

OFX 

' 769 

01401 

OFY 

768 

01400 

PFX 

767 

01377 

PFY 

766 

01376 

RDX 

765 

01375 

RDY 

764 

01374 

RFX 

763 

01373 

RFY 

762 

01372 

XMAX 

761 

01371 

YMAX 

760 

01370 













SYMBOLS AND LOCATIONS FOR SOURCE 

PROGRAM 

FORMAT STATEMENTS 





EFN 

LOC 


EFN 

LOC 


EFN 

LOC 


EFN 

LOC 


EFN 

LOC 

8 ) FK 

500 

01353 

8 ) FM 

502 

01344 

8 ) FU 

510 

01340 

8 )G8 

520 

01336 

8 ) GI 

530 

01332 

8 ) GS 

540 

01325 
















LOCATIONS 

FOR OTHER 

SYMBOLS 

NOT APPEARING 

IN SOURCE 

PROGRAM 





DEC 

OCT 


DEC 

OCT 


DEC 

OCT 


DEC 

OCT 


DEC 

OCT 

1) 

751 

01357 

2) 

687 

01257 

3) 

700 

01274 

6) 

714 

01312 

9) 

748 

01354 

b) 

696 

01270 

C ) Gi 

754 

01362 

C)G2 

755 

01363 

C)G3 

756 

01364 

C ) G4 

757 

01365 

C)G5 

758 

01366 

C)G6 
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PLOTXY - DICTIONARY 

Superscripts are statement numbers in the FORTRAN listing (p. 9) . 

The notation z_ means option z is being used. 

A1 The plotting character. It is either supplied 144 in P(2) by the pro- 
grammer if 1 or specif led 10 ®* 148 to be an asterisk. It Is -written out 251 
in the INFO record. 

BK The third argument in the call list of this subroutine. Its mode (float- 
ing point) differs from the mode the programmer uses (fixed point) in his 
CALL statement, because it is used as a Boolean variable. It is re- 
named 1 ^ BKK to avoid excessive address modification. 

BKEC The Boolean variable name for the number KODE (the sum of the option num- 
bers) . It is "equivalenced" to KK. Branching on single bit positions of 
BKK, each of which represents a different option number being used, 
occurs . 140 > 150 > 160 > 17 0, 180, 190, 200 

C6 An intermediate answer in the calculation of KX (the line count) from 
x^ 320 , an £ jqt (the print position count) from y^? 30 

C7 C6 rounded and truncated to an integer. 322, 332 

DX The scaling parameter for the x scale-factor. It is either supplied 1 ®® in 
P(8) by the programmer if 16, or calculated 121 in PISTUG. 

DY The scaling parameter for the y scale-factor. It is either supplied 1 ®® in 
P(ll) by the programmer if 52 , or calculated 12 ® in PISTUG. 

F10X The value of 10^^^. It is computed 834 " and used 3 ^£ to prepare x^ for 
printing at the left of the first y grid-line if 64 » 

F10Y The value of 1 C)NPWRY # it Is computed® 3 ® and used 344 to prepare y^ for 
printing at the left of the first y grid-line if 64 . 

FX The scaling parameter for the first x grid-label. It Is either calcu- 
lated 1 ® 1 in PISTUG or supplied 1 ® 3 in P(7) by the programmer if 16 . 

FY The scaling parameter for the first y grid-label. It is either calcu- 
lated 1 ®® in PISTUG or supplied 1 ®® in P(l0) by the programmer if 52 . 

I The index of the loop 31 ® -4 ® 8 that processes all pairs of coordinates. 

I is compared 3 ®® to N to provide a branch for the condition that L is not 
a factor of N. 


PLOTXY 


15 



J The counter for the number of blocks ready to be written. It is 

initialized 3004 * 4 - to 1, stepped 400 by 1 after each pair of coordinates is 
processed, tested 330 against L (the maximum number of points in a DATA 
record), and reset 330 to zero when it is equal to L. It replaces 300 the 
value of L if the last DATA record is to be less than L blocks long. J is 
also used as the index of the loop within each output statement 410 ^ 400 j 
430,440 list . gee appendix B, IV. -E. 

K A switch variable set 4 - 00 to 3 and sent- 1 - 04 to PISTUG when scaling param- 

eters for the X array are to be found; set 100 to 2 and sent 108 to PISTUG 
when scaling parameters for the Y array are to be found. 

KBR A manipulation 330 of KSW to create a two-way branch used 340 to control 

picking up Xj_ and Yp for printing in the left-side field if 64 . 

KBR2 A manipulation 340 ^! KSW to create a two-way branch used 347 to control 

picking up Pp if 8. 

KBRM A manipulation 400 of KBW to create a four -way switch variable used 400 to 

control branching to the correct one of four possible output state- 

ments. 410, 420, 430, 440 

KDX A fixed-point name for DX 400 . It is tested- 1 - 04 against zero to recognize 
an error message from PISTUG. 

KDY A fixed-point name for DY . It is tested against zero to recognize 
an error message from PISTUG. 

KHARX An integer representing the range of the X array. It is always calcu- 
lated, either by PL0TXY 18 7 +4r if 16 or by PISTUG^- 2 ^. It is used 23 ® to 
compute KPWRX if 64 . 

KHARY An integer representing the range of the Y array. It is always calcu- 
lated, either by PL0TXY1® 8+ 4 or by PISTUG^- 28 tf 52 . It is used 232 to 
compute KPWRY if 64 . 

KHDG A switch variable initialized^® 3 to 2, and set 2 ® 8 to 1 if 64_. It is 

used 22 ® to control branching to write a special heading over the left- 
side field print-out of the coordinates, 

KK The fixed-point name for K0DE (the sum of the option numbers). It is 
"equivalenced" to BKK and tested!®^ for zero. 

KPWRX The value of 3 - KHARX. If 64, it is computed, 23 ® used 234 to compute 

FIOX, and written out 238 in a special heading. 

KPWRY The value of 3 - KHARY. If 64, it is computed, 232 used 238 to compute 

FIOY, and written out 238 in a special heading. 
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KSW A basic switch variable. It is initialized 1 ^ to zero, set 1 ?® to 1 if 
8, and increased^® by 2 if 64; hence, it may have the value 0, 1, 2, 
or 3. 

KBR » KSW/2 + 1 is computed 33 ® and used 34 ® to control branching to pre- 
pare and Yj_ for print out if 64 . 

KBR2 = KSW + 1 is computed 34 ®_and used 347 to control branching to pick up 
an element of the P array if 8, 

KBRM = KSW + 1 is computed 4 ®® and used 4 ® 2 to control branching to the 
proper output statement. 

KSX One of three scaling parameters for the x scale. It is either calcu- 
lated-^*^- in PISTUG or supplied 4 ® 4 in P(6) by the programmer if 16 . It is 
used 32 ® to compute C6, an intermediate step in the calculation of KX. It 
is written out 231 in the INFO record- It is the number of places the 
decimal point must be moved to the right in the values of OIX and ODX by 
the 1401 subroutine PLOT. 


KSY One of three scaling parameters for the y scale. It is either calcu- 
lated^® in PISTUG or supplied 4 ® 4 in P(6) by the programmer if 52 . It is 
used^SO to compute C6, an intermediate step in the calculation of KY. It 
is written out 234 i n the INFO record. It is the number of places the 
decimal point must be moved to the right in the values of OFY and ODY by 
the 1401 subroutine PLOT. 

KX The array in which the values of KX (number of lines down) are stored 324 
until they are written out in a DATA record by one of four output state- 
ments • ^10j 420 ^ 430^ 440 


KY 


L 


M 


N 


The array in which the values of KY (number of print positions to the 
right) are stored 334 until they are written out in a DATA record by one of 
four output statements. 410,420^,430^440 


The maximum number of blocks possible in one DATA record. It is com- 
puted 3 ^® and used 33 ® to test J and to write out 4 !®* 42 ®* 43 ®* 44 ® DATA 
record. If there are less than the maximum number of blocks to be written 
in the last DATA record, L is set 3 ® 2 equal to J before the write-out. 

The frequency of the x grid-lines. It is either supplied 434 in P(3) by 
the programmer if 2_ or set!® 7 * I* 3 ® equal to 10. It is written out 234 in 
the INFO record. 


The number of points to be plotted. It is supg 
in P(l).^dused in both calls to PISTUG. 424 * 42 
of executions of the loop 3 !®*" 4 ®® that processes 
the DATA records. It is tested 3 ®® against I to 
condition that L is not a factor of N. 


lied!®® by the programmer 
® It is the total number 
all points and writes out 
cause branching on the 
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HF1 


m i 

OIK 

ODY 

OFX 

OFY 

P 

PFX 

PPY 

POW 

RDX 

RDY 

REX 


The length of the first field in one block of a DATA record. The NF1 
characters in a block are to be printed to the left of the first y grid- 
line on the same line with the point whose KX and KY values are in the 
block. If NF1 = 0, blanks are_to be printed. HF1 is initialized 1 ®^/ 17 ® 
to zero, set 174 equal to 6 if 8, and Increased 2 ® 4 by 13 if 64 . It is 
used®®® to compute L, the maximum number of blocks per record, and is 
written 2 ® 1 out in the INFO record. 

The frequency of the y grid-lines. It is_either set 1 ®®/ 1 ®® to 10 or 
supplied^® 4 in P(4) by the programmer if 4, and it is written out 2 ® 1 
in the INFO record. 


A number written out 2 ® 1 representing the x scale-factor. It is com- 
puted 2 ®® by multiplying the scaling parameter DX (always an integer) by 
10“® to put it in the form required for the INFO record. 

A number written out 2 ® 1 representing the y scale-factor. It is com- 
puted 2 ®®' 1 ' 2 by multiplying the scaling parameter DY (always an integer) 
by 10"® to put it in the form required for the INFO record. 


A number written out 2 ® 1 representing x-scale starting -value. It is com- 
puted®®®'*" 1 by multiplying the scaling parameter FX (always an integer) by 
10“® to put it in the form required for the INFO record. 


A number written out 2 ® 1 representing the y-scale starting-value- It is 
computed 2 ®® - *"® by multiplying the scaling parameter FY (always an integer) 
by 10~® to put it i,n the form required for the INFO record. 


An array name, the fourth argument in the call list of this subroutine. 
The number of points to be plotted must be the first element of this 
array and the other elements are any data required by the particular op- 
tions being used. (See ref. 1, p,7) 

An intermediate value found 1 ® 7 "*" 2 in the calculation of KHARX. 

An intermediate value found 1 ®®'*' 2 In the calcination of KHARY. 

The array that has as its i 1 * 1 element the six BCD characters_associated 
with Xj_. It is supplied® 4 ® in P(l+ll) by the programmer if £ and written 
out in the DATA record. 42 ®/ 44 ® 

The true value of the x scale-factor. It is computed and used 1 ® 1 el seq. 
to compute KHARX if 16 . 

The true value of the y scale-factor. It is computed and used 1 ®® el seq. 
to compute KHARY if 32 . 

The true value of the x-scale starting -value. It is computed and 
used 1 ® 7 e ’ t: ' se l* to compute KHARX if 16. 
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RFY 

X 

XMAX 

XOW 

Y 

YMAX 

YOW 


■ I Ml I ■■■ II IMI 


I II II III III II I I III 


iiiii in mi i i i iiiii i i ii i mi ii 


The true value of the y-scale star ting -value. It is computed and 
US ed.l99 et seq. to compute KHARY if 32. 

An array name, the first argument in the call list of this subroutine, 
containing the values of the variable to be plotted on the vertical 
scale. It is sent 121 to PISTUG for scaling, unless 16. The elements 
must be in sequence, either increasing or decreasing. 

An intermediate value found^- 8 ^ sec l* during the calculation of KHARX. 

The array in which the transformed values of X± are held^2 until they 
are written out^30,440 the DATA record if 64 . 

An array name, the second argument in the call list of this subroutine, 
containing the values of the variable to be plotted on the horizontal 
scale. It is sent 128 to PISTUG for scaling, unless 52 . 

An Intermediate value found^- 88 e ^ sec l* during the calculation of KHARY. 

The array in which the transformed values of Yi are held^^ until they are 
written out^0,440 i n the DATA record if 64. 
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PLOTMY - (FORTRAN II) 


Although more than one curve can be plotted on the same grid using FLOTXY, 
they will all use the same plotting character- When a different plotting char- 
acter for each curve is desired, the programmer uses PLOTMY. The CALL statement 
lists the names of the arrays to be plotted; an array containing the number KODE, 
the number of curves, and possibly other information - depending upon the options 
(PLOTXY I.-D.) used and the Variation (PLOTMY I.-C.) chosen; and finally an array 
whose first element indicates the Variation chosen and whose other elements are 
any additional data required for the options being used. 

PLOTMY may call an auxiliary subroutine PISTUG to compute scaling param- 
eters, then it writes an INFO record specifying the appearance of the grid to 
the 1401 subroutine PLOT. This is followed by one or more packed DATA records 
that specify the plotting character as well as the location of each point in 
terms of KX (the number of lines down from the first x grid-line) and KY (the 
number of print positions to the right of the first y grid-line). 

I. Input to this routine is supplied via the four arguments of the CALL, as 
follows : 


CALL PLOTMY (XDOWN, YACROS, KECK, P) 

A. XDOWN and YACROS are floating-point variable array names. In general, 
the values in both arrays must be between -10® and 10®. (See appen- 
dix B,VI. , for precise limitations.) The subroutine SCALE is supplied 
to handle arrays whose element size is not known. (SCALE should be 
called before calling PLOTMY.) There is no restriction on the order of 
the elements; however, the elements of the array XDOWN .are destroyed. 

The minimum DIMENSIONS of these arrays must be as shown in reference 1, 
page 11. 

B. KKK is a fixed-point variable array name, whose minimum DIMENSIONS 
must be either 14, or two more than twice the number of curves to be 
plotted, whichever is larger. The elements in this array are determined 
by the choice of options (i.-E.) and Variation (see I.-C.) made by the 
programmer. 

KKK(l) is the KODE. Seven different options are offered to the pro- 
grammer. Each has a number associated with it. The sum of the option 
numbers Is the KODE, and 0 < KODE < 127 . 

KKK(2) = KN, the number of curves to be plotted. If KN > 6, option 1 
must be used. 

KKK(3) is supplied as specified by the Variation selected by the pro- 
grammer . 
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KECEC(4) , KKK(6) , « • . are supplied only if option 1 is used. 

KKK(5 ) j KKK(7), . . . are supplied only if Variation III is selected. 
G. Variations 

Frequently a set of curves to be plotted is defined by several sets 
of ordinates associated -with the same set of abscissas (each value of an 
independent variable determining one value of each of several dependent 
variables). Variation I is for such multiple curves, with the single set 
of abscissas placed in the XDOWN array and the multiple sets of ordinates 
in YACROS. For complete flexibility, Variation II permits the inter- 
change of axes; that is, the single set of abscissas in YACROS, the ordi- 
nates in XDOWN. In these tvo Variations each curve has the same number 
of points (MFCS), and it is supplied to the subroutine in KECK(3)* 

Variation III is for general multiple curves, in which each curve 
may have any number of points . The number of points for the first curve 
is supplied in KKK(3), for the second curve in KKK(5); in general, the 
number of points for the curve must be in KKK(2*J+l) . 

D. P is a floating-point variable array name: 

P(l) = 1. for Variation I. 

P(l) = 3. for Variation II. 

P(l) = 5. for Variation III. 

The DIMENSION and remaining elements of the P array are determined 
by the particular options being used (see ref. 1, p. ll) . 
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E. Options. The sample plots in appendix A illustrate the effects of 
using the various options . 


TIMBER 


OPTION 


1 


2 


Choice of plotting characters . - The plotting characters 
+, 0 , x, =; C will he used for the first six curves un- 
less the programmer uses this option. When he does^ he 
must supply the desired plotting characters (which may he 
any FORTRAN character except the minus sign) in the KKK 
array. The octal representation of the plotting character 
for the first curve must he in the leftmost 6 hits of 
KKK(4) . for the second curve in KKK(6) and; in general; for 
the curve in KKK(2*J+2) . If more than six curves are 
to he plotted^ the programmer must use option 1; and must 
supply all the plotting characters. 


Choice of x grid-line frequency . - An x grid-line will he 
printed every 10 line-spaces helow the first one unless the 
programmer uses this option. When he does^ the desired fre- 
quency (M) of x grid-lines must he supplied to PLOTMI in 
P(3) . If M a 0; only the first and final x grid-lines will 
he printed. 


4 Choice of y grid-line frequency . -Ay grid-line will he 

printed every 10 print positions to the right of the first 
unless this option is used. If it is, the desired fre- 
quency (N) of y grid-lines must he supplied in P(4), If 
N ss 0; only the first (leftmost) y grid-line will he 
printed. 


16 Choice of x scale . - PLOTMT will supply the three scaling 

parameters that specify the x-scale starting -value; the 
scale-factor (value of one line-space); and the x grid- 
lahels unless this option is used. If it is used; the pro- 
grammer must specify his own scale by supplying the three 
parameters KSX; FX; and DX in P(6); P(7); and P(8) , re- 
spectively. These parameters are defined hy: 


FX is the desired starting -value for the x scale 
multiplied hy 10 U . 


DX is the desired scale-factor for one line-space 
multiplied by 10^. DX must he positive. 

U must he chosen so that FX and DX are integers and 
0 < U < 6. 


KSX is equal to 6 - U. 
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NUMBER OPTION 

32 Choice of y scale , - FLOTMY will supply the three scaling 

parameters that specify the y-scale star ting -value, the 
scale-factor (value of one print position) and the y grid- 
labels unless the programmer uses this option. If he does, 
he must specify his own scale by supplying the three param- 
eters, KSY a FY, and DY in P(9), P(lO), and P(ll), respec- 
tively. These are defined by: 

FY is the desired starting -value multiplied by 10^. 

DY is the desired scale-factor for one print position 

multiplied by 10*^. 

V must be chosen so that FY and DY are integers and 

0 < V < 6. 

KSY is equal to 6 - V. 

64 C hoice of printing x-coordinate at the left of the plot . - 

The 29 print positions to the left of the first y grid-line 
are usually blank except for x grid-labels every tenth line 
in positions 20 to 28. If option 64 is used, the x-ordinate 
of each point is printed immediately to the left"*- of the 
grid-labels on the same line as the point. 

8 Choice of printing message at left of plot . - The 29 print 

positions to the left of the first y grid-line are usually 
blank except for x grid-labels every tenth line in posi- 
tions 20 to 28. If option 8 is used, the 6 BCD characters 
in P(l2) are printed immediately to the left-** of the x grid- 
labels on the same line on which the first point is printed, 
the six characters in P(l3) are printed with the second 
point, etc. Specifically, the 6 characters in P( I+ll) are 
printed on the same line on which the I^ value of the XDOWN 
array is plotted. 


Ilf option 8 and 64 are used together, the BCD characters will print to the 
left of the ordinate. 
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II. Output from this routine consists of one INFO record and one or more 
packed DATA records. Sample tape records appear in appendix A. As a result 
of the execution of this routine, the elements of the array XDOWN are de- 
stroyed. 

A. INFO record. 

The contents of the INFO record are the information needed by the 
1401 subroutine PLOT to determine the construction of the grid, the con- 
struction of the grid-labels, and the format of the DATA records that 
follow. Specifically, the contents of the INFO record consist of the two 
characters, PI, followed by 10 fields, whose contents are: 

(Field width is shown in parentheses.) 

a. (l) This character is always a minus sign, and signifies that the 
plotting characters to be associated with each point will be found in 
the DATA records. 

b. (l) This digit is the decimal -point-shift code, KSX, for OFX and 

ODX. (items d. and f. in this list.) 

c. (l) This digit is the decimal -point-shift code, KSY, for OFY and 

ODY. (items e. and g. in this list.) 

d. (9) The number (OFX) in this field is in the form iO.xxxxxx and 
represents the x-scale star ting -value. 

e. (9) The number (OFY) in this field is in the form iO.xxxxxx and 
represents the y-scale star ting -value. 

f. (9) The number (ODX) in this field is in the form IO.xxxxxx and 
represents the x scale-factor (value of one line-space) . 

g. (9) The number (ODY) in this field is in the form iO.xxxxxx and 
represents the y scale-factor (value of one print position) . 

h. (3) This integer is the frequency, M, of the x grid-lines. 

i. (3) This integer is the frequency, N, of the y grid-lines. 

j. (3) This integer specifies the length, NF, of the first field in 
each block of information in the following DATA records. NF may be 
zero. 
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B. DATA records. 


The contents of a DATA record (after the first two characters PD) 
are successive blocks of information needed by the 1401 subroutine PLOT 
to position one or more plotting characters in the same line and print 
the left-side field of that line. Within one record, only complete 
blocks are permitted, and each block is terminated by an end-of -block 
character or delimiter - the letter E. One record may hold as many as 
16 blocks. Within one block, the first W characters are those to be 
printed in the left-side field (to the left of the first y grid-line) un- 
less W =* 0; in this case, blanks are to be generated by the 1401. The 
next three characters are KX, the number of line-spaces below the first 
x grid-line at which all the points in that block are to be printed. The 
next three characters are a value of KY, the number of print positions 
to the right of the first y grid-line, in which a point is to be printed. 
The plotting character itself follows KY. Since additional points may 
fall on the same line, further values of KY, each with its associated 
plotting character, may follow the first. The plotting character for 
the last value of KY in the block is followed by the delimiter E. 

III. After initializing switches, the routine inspects the input variable 
KODE, bit by binary bit, branching on each of the seven bit positions to 
set necessary output switches, and either supplies standard values or picks 
up programmer -provided specifications from the KKK and P arrays. It 
branches also on Variation number for additional initializing. It may call 
an auxiliary routine PISTUG (p. 44) for scaling either or both coordinate 
arrays. The INFO record (H.-A.) is written, and the routine constructs 
the output FORMAT statement for the packed DATA records. 

The XDOWN array is searched for a minimum value and the corresponding 
value of KX (number of lines down) is computed and stored in an output 
buffer area. KY*s (the number of print positions across) are computed from 
corresponding YACROS values; each is sent to the buffer with its associated 
plotting character; the last is followed by the delimiter. The minimum 
element of the XDOWN array is then overwritten with a large number, and a 
new search for a minimum is made. Whenever the output buffer area is 
filled, a DATA record is written. This process continues until all points 
have been processed; control is then returned to the calling program. 
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DO 

a 


SUBROUTINE PLOTMYt X, Y,K,P ) 


C 


C 

C 


C 

B 

B 

B 

B 


C 


c 

c 

B 


C 

C 

B 

8 


SUBROUTINE PLOTMY { X, Y» K, P ) 

EQUIVALENCE ( BKK.KKJ , ( A,KPC| , CFUSEt 5) tFUSEMI » i B,KB) 

DIMENSION X(l)«Ytl)«K(l),P(l),A(6) ,KPC( 6 ) , FUSE! 21 ) , FUSEMt 16) , B t 64 ) 

1» KB ( 64 ) 

RENAME INPUT 

100 KK-K(l) 

101 KN=Kt2) 

102 NPTS*K ( 3) 

SET UP FORMAT STATEMENT FOR DATA OUTPUT 

103 FUSE l 1 ) =740230472473 

104 FUSE! 2) * 020074606060 

105 FUSE ( 3) = 606060606060 

106 FUSE t 4 ) * 603103736060 

107 DO 108 1=1*16 

108 FUSEMll) = FUSEl 3) 

109 FUSE( 21 ) = 6H Al)) 

INITIALIZE SWITCHES, COUNTERS* CONSTANTS 

110 KEX IT*i 

111 KTL*1 

112 KSW64* 1 

113 KSW8-1 

114 KSW1*1 

116 KSW3=1 

117 KWAOO*2 

118 KLA0D*4 

119 NF1*0 
M = 10 
NY=10 

OPTION 1 1 GET STD. P.C.S OR GO TO 150 

120 IF(BKK*1000000) 122,122,150 
1 22 All)* 1H* 

A( 2 ) 3 1H+ 

At 3 ) = 1H0 
A ( 4 ) 3 1HX 
At 5 ) 3 1H- 
A t 6 ) = 1H0 
K 1 4 ) * KPCtl) 

K t 6 ) = KPC 1 2 ) 

K t 8 ) * KPC ( 3 ) 

K 1 10 ) *KPC ( 4 ) 

K( 12) =KPC { 5 ) 

Kt 14) *KPC t 6 ) 

GO TO 200 IF KODE IS ZERO 

140 IF(KK) 150,200,150 

INITIALIZE FOR OPTION 64 OR GO TO 160 

150 IF(8KK*100000000> 160,160,152 
152 KSW64=2 

154 FUSEt 3 ) =602607330073 
156 NF 1=7 

158 KLADD*KLADD+7 

159 KWADD»KWADD+1 


C 

C 



PLOTMY 




SUBROUTINE PLOTMY( X, Y,K,P) 

c 

B 160 IF(6KK*10000000) 170,170,162 
162 KSW8=2 

B 164 FUSE(2)=020074210673 
166 NFl=NFi+6 

168 KLADD=KLA0U+6 

169 KHADU=KWADD+1 
C 

8 170 IF( BKK*2000000) 180,180,172 
172 M=P(3)+SIGNF(.5,P(3) ) 

C 

B 180 I F ( BKM4000000) 200,200,182 
182 NY=P(4)+SIGNF( • 5 , P ( 4 ) ) 

C 

200 IF ( P ( 1 )-2 • 5 ) 202,210,210 
202 KTL=KN 
204 KSH1=2 
208 GO TQ 228 
C 

210 I F { P ( 1 ) -4 • ) 220,220,214 
214 NPTbT-0 
216 DO 218 1=1, KN 

218 NPTST=NPTST+K(2*I+1) 

219 GO TO 230 
C 

220 KTIMES=KN-1 

221 DO 225 1 = 1 » KT IMES 

222 MM=NP TS* I 
1222 K ( 2* I + 3 ) = NPTS 

223 DO 225 11=1, NPTS 

224 L = MM+ 1 1 

225 Y(L) = YUI> 

226 KSW3=2 
C 

228 NPTST=KN*NPTS 
C 

B 230 I F ( BKK*20000000)250,250,23l 

231 KSX=P(b)+SlGNF( .5,P(6) ) 

232 FX = P ( 7 J 

233 DX=P ( 8 ) 

234 RFX=FX*(10.**(KSX-6) ) 

235 R0X=DX*(10.**(KSX-6 ) > 

236 PFX=RFX+100.*RDX 

237 XMAX=MAX1F (ABSF(RFX) ,ABSF(PFX) ) 

238 KHARX=X I NTF ( LQGF ( XMAXJ/2. 302585+40 
C 

B 240 I F ( BKK*400Q0000) 270,270,242 
242 KSY=P(9)+SIGNFl.5,P(9) ) 

244 F Y=P ( 10 ) 

246 D Y=P ( 11) 

248 GO TO 300 
C 

250 GO TO (256,252) ,KSW1 
252 NS X=NP T S 
254 GO TO 260 
256 NSX=NP I ST 


INITIALIZE FOR OPTION 8 OR GO TO 170 

IF OPTION 2, USE P ( 3 ) 

IF OPTION 4, USE P(4) 

INITIALIZE FOR V.I OR GU TO 210 

INITIALIZE FOR V.III OR GO TO 22 0 
INITIALIZE FOR V.II 


FOR V.I AND V.II 
IF OPTION 16, USE P ( 6-8 ) 


)-40 

IF OPTION 32, USE P(9-ll) 
NSX DIFFERS FOR V. I 



PLOTMY 


t\3 

CD 


SUBROUTINE PLOTMY ( X, Y,K, P ) 


C SCALE X AND TO 240 OR ERROR EXIT 

260 CALL P ISTUG ( 1 , NSX , X f OX, FX , KSX, KHARX I 
262 KOX=DX 

264 IF(KDX) 240,900,240 

C NSY DIFFERS FOR V. 1 1 

270 GO TO (276,272) ,KSW3 
272 NS Y=NP T S 
274 GO TO 280 
276 NSY=NPTST 

C SCALE Y AND TO 300 OR ERROR EXIT 

280 CALL PISTUG l 2, NSY, Y, DY , FY ,KSY ,KHARY ) 

282 KQY=DY 

284 IF1KDY) 300,900,300 

C WRITE 0PT64 HEADING OR GO TO 310 

300 GO TO 1310, 302), KSW64 
302 KPWRX=3-KHARX 
304 F10X=10.**KPWRX 
306 WRITE OUTPUT TAPE 6,500,KPWRX 
500 FOR MAT(2 HPT/2 HPT/2 HP T,15X,3HX*E, 12) 

C WRITE THE INFO RECORD 

310 RO* 5.6-7 

311 SFT=1 • 6-6 

312 OOX=DX#SFT+SIGNF1RO,DX) 

313 OFX=FX*SFT+SIGNFlRD,FX) 

314 ODY=DY*SFT+SIGNF(RD,DY) 

315 0FY=FY*SFT+S1GNFIRD,FY> 

316 XMF=I 10.**( 6-KSX) )/DX 

317 YMF=( 1U.**(6-KSY) )/DY 

318 WRITE OUTPUT TAPE 6,502,KSX,KSY,0FX,OFY,0DX,0DY,M,NY,NFl 
C 

502 FORMAT 13HP1-,2I1,4F9.6,3I3) 

C 

C INITIALIZE FOR DATA OUTPUT 

C 

C FINISH FORMAT STATEMENT 

320 DO 322 I = 1, KTL 
B 322 FUSEM(I) * 310373210173 

C COMPUTE BUFFER CONSTANTS AND SET CTRS 

330 KLC=KLADD+4«KTL 
332 NBR= 1 30/KLC 
334 KW =KWAD0+2*KTL 

336 NBW=KW*NBR 

337 J = 1 

338 KB W= 1 

339 KBBC* l 

340 KLAST=2*KN+1 

C ILIM DIFFERS FOR V.I 

342 GO TO 1348, 344), KSW1 
344 ILIM=NPTS 
346 GO TO 350 
348 ILIM*NPTST 

C READY TO FIND SMALLEST X 

350 XM IN* 1 . E15 
352 1MIN=1 
C 
C 



PLOTMY 


SUBROUTINE PLOTMY ( X » Y » K , P ) 


CO 

tO 


C 

360 DO 368 1 = 1, ILIM 
362 IF(XMIN-XU)) 368,368,364 
364 XM1N=XU) 

366 I M I N= I 
368 CONTINUE 
C 

370 GO TU (380,372) ,KSW8 
372 B t KB W ) =P ( IMIN+11) 

374 KB W=KB h + 1 
C 

380 GO TU (390,382) ,KSW64 
382 B(KBW)=F10X*XMIN 
384 KBW=KBW+1 
C 

390 C6 = XM I N* XMF 
392 C7 =1NTF(C6+SIGNF( .5,C6) ) 

394 KBt K6W)=( DX*C7-FX) /DX+.5 
396 KBW=KB W+ 1 
C 

398 Xi IMIN) = 1.E15 
C 
C 
C 

400 DO 446 1M= 1 , KTL 
402 LL=IMlN+( IM-1)*NPTS 
404 YS=Y(LL) 

406 C6 = YS* YMF 

408 C7=INTF(C6+SIGNF( .5,C6) ) 

410 KB(K8W)= (DY*C7-FY)/DY+.5 
412 KBW=K6h+1 
C 

420 GO TO 1422, 440 ) , KSW 1 
422 I K = 0 

424 DO 4301L=3,KLAST,2 
426 I K= IK+K ( IL) 

428 IF l IK- IM IN ) 430,436,436 
430 CONTINUE 
C 

432 WRITE OUTPUT TAPE 6,600, IMIN 
600 FORMAT ( 2HPL , 42H NPTS IS MISSING 
434 GO TO 900 
C 

436 KPC I * K ( IL+1) 

438 GO TO 442 
C 

440 KPC I = K ( 2* IM+2 ) 

C 

442 KB ( KB w ) =K PC I 
444 KB W = K8 w + 1 
446 CONTINUE 
C 

448 B ( KBW ) = IHE 

449 KB W=KB W+ 1 
C 

450 IF(NBW-KBW) 460,452,452 


FIND SMALLEST X AND ITS SUBSCRIPT IMIN 

IF OPTION 8, MOVE PUMlN+llI TO BUFFER 
IF OPTION 64, MOVE X ( IMIN) TO BUFFER 
COMPUTE KX AND MOVE TO BUFFER 

OVERWRITE SMALLEST X 

ASSOCIATED KYS AND P.C.S GO TO BUFFER 

IF NOT V. I, SEARCH FOR P.C. 

IF NOT FOUND, ERROR EXIT 
OR WRONG IN KARRAY.IMIN= ,13) 

P.C. FOUND FOR V.II AND V.III 

P.C. FOR V. I 

MOVE PLOT CHAR. TO BUFFER 

DELIMITER TO BUFFER 

GO TO 460 IF BUFFER FILLED 



PLOTMY 


SUBROUTINE PLOTMY { X, Y, K , p ) 


O 


C 

452 

C 

454 

456 

458 

C 

460 

462 

C 

464 

466 

C 

468 

C 

470 
C 

471 

472 
C 

900 


I F ( IL IM-J ) 470,470,454 


GO TO 470 IF ALL POINTS COMPUTED 
STEP COUNTERS AND GO BACK TO 350 


J=J + 1 

KB8C=KB8C+1 
GO TO 350 

WRITE A DATA RECORD 

KBL-KW»KBBC 

WRITE UUTPUT TAPE 6, FUSE, ( B ( N ) , N* 1 , KBL ) 

RESET BUFFER COUNTERS 


KBW=1 

KBBC=0 

GO TO 1 452 , 900 ) , KEXIT 

IF(KBBC) 471,900,471 

KEX IT=2 
GO TO 460 

RETURN 

END (1,0, 0,0 ,0,0, 0,0, 0,1, 0,0, 0,0,0) 


GO TO 452 OR QUIT 

IF ALL RECORDS WRITTEN, QUIT 

IF NOT f SET EXIT TO QUIT AND GO TO 460 



PLOTMY 


SUBROUTINE PLOTMY ( X , Y, K , P ) 


STORAGE NOT USED BY PROGRAM 

DEC OCT DEC OCT 

1125 02145 32561 77461 


STORAGE LOCATIONS FOR VARIABLES APPEARING IN DIMENSION AND EQUIVALENCE STATEMENTS 



DEC 

OCT 


DEC 

OCT 


DEC 

OCT 


DEC 

OCT 


DEC 

OCT 

A 

1123 

02143 

BKK 

1124 

02144 

B 

1096 

02110 

FUSEM 

1113 

02131 

FUSE 

1117 

02135 

KB 

1096 

02110 

KK 

1124 

02144 

KPC 

1123 

02143 








STORAGE LOCATIONS FOR 

VARIABLES NOT APPEARING 

IN COMMON, DIMENSION, OR EQUIVALENCE 

STATEMENT 




DEC 

OCT 


DEC 

OCT 


DEC 

OCT 


DEC 

OCT 


DEC 

OCT 

C6 

1032 

02010 

C7 

1031 

02007 

DX 

1030 

02006 

DY 

1029 

02005 

F10X 

1028 

02004 

FX 

1027 

02003 

FY 

1026 

02002 

11 

1025 

02001 

IK 

1024 

02000 

IL I M 

1023 

01777 

IL 

1022 

01776 

IMIN 

1021 

01775 

IM 

1020 

01774 

I 

1019 

01773 

J 

1018 

01772 

KBBC 

1017 

01771 

KBL 

1016 

01770 

KBW 

1015 

01767 

KUX 

1014 

01766 

KDY 

1013 

01765 

KEXIT 

1012 

01764 

KHARX 

1011 

01763 

KHARY 

1010 

01762 

KLADD 

1009 

01761 

KLAST 

1008 

01760 

KLC 

1007 

01757 

KN 

1006 

01756 

KPC I 

1005 

01755 

KPWRX 

1004 

01754 

KSW1 

1003 

01753 

KSW3 

1002 

01752 

KSW64 

1001 

01751 

KSW8 

1000 

01750 

KSX 

999 

01747 

KSY 

998 

01746 

KTIMES 

997 

01745 

KTL 

996 

01744 

KWADD 

995 

01743 

KW 

994 

01742 

LL 

993 

01741 

L 

992 

01740 

MM 

991 

01737 

M 

990 

01736 

NBR 

989 

01735 

NBW 

988 

01734 

NF1 

987 

01733 

NPTS 

986 

01732 

NPTST 

985 

01731 

NS X 

984 

01730 

NSY 

983 

01727 

NY 

982 

01726 

ODX 

981 

01725 

ODY 

980 

01724 

OFX 

979 

01723 

OFY 

978 

01722 

PFX 

977 

01721 

RD 

976 

01720 

RDX 

975 

01717 

RFX 

974 

01716 

SFT 

973 

01715 

XMAX 

972 

01714 

XMF 

971 

01713 

XMIN 

970 

01712 

YMF 

969 

01711 

YS 

968 

01710 




SYMBOLS AND LOCATIONS FOR SOURCE 

PROGRAM 

FORMAT STATEMENTS 





EFN 

LOC 


EFN 

LOC 


EFN 

LOC 


EFN 

LOC 


EFN 

LOC 

8 ) FK 

500 

01657 

8 ) FM 

502 

01652 

8) 10 

600 

01646 










LOCATIONS 

FOR OTHER 

SYMBOLS NOT APPEARING 

IN SOURCE 

PROGRAM 





DEC 

OCT 


DEC 

OCT 


DEC 

OCT 


DEC 

OCT 


DEC 

OCT 

1) 

951 

01667 

2) 

869 

01545 

3) 

895 

01577 

6) 

919 

01627 

9) 

944 

01660 

B) 

879 

01557 

C)G2 

954 

01672 

C)G3 

955 

01673 

C)G5 

956 

01674 

C)G6 

957 

01675 

C )G7 

958 

01676 

C )G8 

959 

01677 

C)G9 

960 

01700 

UGA 

961 

01701 

C)GB 

962 

01702 

C ) GC 

963 

01703 

C)GD 

964 

01704 

C ) GF 

965 

01705 

0200 

966 

01706 

0201 

967 

01707 

D) 10J 

280 

00430 

0J11H 

653 

01215 

D) 1 IS 

784 

01420 

D 1 1 1 T 

788 

01424 

D) 11V 

814 

01456 

D) 121 

828 

01474 

D) 126 

861 

01535 

D ) 206 

186 

00272 

D) 20C 

239 

00357 

D) 20 J 

282 

00432 

D) 20N 

316 

00474 

D) 200 

322 

00502 

D ) 2 IK 

706 

01302 

D)30N 

315 

00473 

0)300 

321 

00501 

0)408 

207 

00317 

D) 410 

444 

00674 

D ) 526 

860 

01534 

E ) 1 H 

655 

01217 

E ) 1 J 

668 

01234 

E)1L 

746 

01352 

E ) 20 

820 

01464 

E ) 326 

863 

01537 












LOCATIONS 

OF NAMES 

. IN TRANSFER 

VECTOR 







DEC 

OCT 


DEC 

OCT 


DEC 

OCT 


DEC 

OCT 


DEC 

OCT 

EXPI2 

0 

00000 

LOG 

1 

00001 

PISTUG 

2 

00002 

(FIL) 

4 

00004 

I STH) 

3 

00003 


ENTRY POINTS TO SUBROUTINES NOT OUTPUT FROM LIBRARY 


PLOTMY 


ro 


SUBROUTINE PLOTMY { X f Y, K , P ) 


EXP { 2 

LOG 

EXTERNAL 

PISTUG (FIL) 

FORMULA NUMBERS WITH 

( STH) 

CORRESPONDING 

INTERNAL FORMULA 

NUMBERS AND OCTAL 

LOCATIONS 



EFN 

IFN LOC 

EFN 

IFN LOC 

EFN 

IFN LOC 

EFN 

IFN LOC 

EFN 

IFN 

LOC 

100 

8 00123 

101 

9 00125 

102 

10 00127 

103 

11 00131 

10* 

12 

00133 

105 

13 00135 

106 

1* 00137 

107 

15 001*1 

108 

16 001*2 

109 

17 

001*6 

no 

18 00150 

111 

19 0015* 

112 

20 00156 

113 

21 00162 

11* 

22 

00166 

116 

23 00172 

117 

2* 00176 

118 

25 00200 

119 

26 00202 

120 

29 

00210 

122 

30 0021* 

1*0 

*2 002** 

150 

*3 002*6 

152 

** 00253 

15* 

*5 

00257 

156 

*6 00261 

158 

*7 00263 

159 

*8 00266 

160 

*9 00273 

162 

50 

00300 

16<V 

51 0030* 

166 

52 00306 

168 

53 00311 

169 

5* 0031* 

170 

55 

00320 

172 

56 00325 

180 

57 00337 

182 

58 003** 

200 

59 00360 

202 

60 

0036* 

20* 

61 00366 

208 

62 00372 

210 

63 00373 

21* 

6* 00*00 

216 

65 

00*02 

218 

66 00*07 

219 

67 00*1* 

220 

68 00*15 

221 

69 00*20 

222 

70 

00*33 

1222 

71 00*37 

223 

72 00**1 

22* 

73 00**6 

225 

7* 00*53 

226 

75 

00*66 

228 

76 00*75 

230 

77 00503 

231 

78 00510 

232 

79 00522 

233 

80 

0052* 

23* 

81 00526 

235 

82 005*2 

236 

83 00556 

237 

8* 00562 

238 

85 

00575 

2*0 

86 00613 

2*2 

87 00620 

2** 

88 00632 

2*6 

89 0063* 

2*8 

90 

00636 

250 

91 00637 

252 

92 006*2 

25* 

93 006** 

256 

9* 006*5 

260 

95 

006*7 

262 

97 00661 

26* 

98 00670 

270 

99 00675 

272 

100 00700 

27* 

101 

00702 

276 

102 00703 

280 

103 00705 

282 

105 00721 

28* 

106 00730 

300 

107 

00732 

302 

108 00735 

30* 

109 007*0 

306 

110 007*6 

310 

112 0075* 

311 

113 

00756 

3L2 

1L* 00760 

313 

115 00770 

31* 

116 01000 

315 

117 01010 

316 

118 

01020 

317 

119 01032 

318 

120 010** 

320 

' 122 01072 

322 

123 01075 

330 

12* 

01101 

332 

125 01106 

33* 

126 0111* 

336 

127 01121 

337 

128 01125 

338 

129 

01127 

339 

130 01133 

3*0 

131 01135 

3*2 

132 011*2 

3** 

133 011*5 

3*6 

13* 

011*7 

3*8 

135 01150 

350 

136 01152 

352 

137 0115* 

360 

138 01160 

362 

139 

01165 

36* 

1*0 01172 

366 

1*1 0117* 

368 

1*2 01200 

370 

1*3 0120* 

372 

1** 

01207 

37* 

1*5 01211 

380 

1*6 01221 

382 

1*7 0122* 

38* 

1*8 01227 

390 

1*9 

01235 

392 

150 012*0 

39* 

151 012*7 

396 

152 01263 

398 

153 01270 

*00 

15* 

01272 

*02 

155 01303 

*0* 

156 0131* 

*06 

157 01316 

*08 

158 01321 

*10 

159 

01330 

*12 

160 013** 

*20 

161 01350 

*22 

162 01353 

*2* 

163 01355 

*26 

16* 

01361 

*28 

165 0136* 

*30 

166 01377 

*32 

167 01*02 

*3* 

169 01*1* 

*36 

170 

01*15 

*38 

171 01*17 

**0 

172 01*21 

**2 

173 01*25 

*** 

17* 01*27 

**6 

175 

01*33 

**8 

L76 01**2 

**9 

177 01*** 

*50 

178 01*51 

*52 

179 01*57 

*5* 

180 

01*65 

*56 

181 01*70 

*58 

182 01*73 

*60 

183 01*75 

*62 

18* 01501 

*6* 

189 

01516 

*66 

900 

190 01521 
195 015*1 

*68 

191 01523 

*70 

192 01526 

*71 

193 01530 

*72 

19* 

01533 



33 



152-159 162-169 172 182 



PLOTMY 


372 382 



PLOTMY 
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PLOTMY - DICTIONARY 


A 

B 

BKK 

C6 

C7 

DX 

DY 

FlOX 

FUSE 


Superscripts are statement numbers in the FORTRAN listing (p. 26) . 

The notation _z means option z is being used. 

An array that holds the six standard plotting characters. It is equi va- 
lence! to KPC. The array is set up as A and used as KPC122 et seq un- 
less 1. 

The array name of the output buffer area (or buffer) in which DATA 
records are assembled. It is equivalenced to KB. The name B is 
used372>382,448 bo move the left-side field characters from the P array 
if £3, x^ if 64, and the delimiter E (end-of -block character). The B 
array is written out 462 into the DATA record. 


The Boolean variable name for KODE (the sum of the option numbers). It 
is "equivalence!" to KK. Branching on single bit positions, each 'of 
which represents a different option number being used, occurs. 120, 150, 
160,170,180,230,240 


An intermediate answer in the calculation of KX 360 from xq and KY^ 06 
from y-£. 

C6 rounded and truncated392,408 to an integer. 


The scaling parameter for the x scale-factor. It Is either supplied 233 
in P(8) by the programmer if 16_ or calculated 260 in PISTUG. It is used 
to compute XMF 316 and each value of KX? 94 If' 16, it Is used 235 to 
compute REK. See also KDX^ 0DX. 


The' scaling parameter for the y scale-factor. It is either supplied 2 ^ 6 
In P(ll) by the programmer if 52 , or calculated 260 in PISTUG. It is 
used to compute YMF 3 ! 7 and each value of KY^ 10 See also KDY, 0DY. 


The value 30 ! of io^^RX. It Is used to prepare 362 X^ for printing at 
the left of the first y grid-line if 64. 


The array (DIMENSION is 2l) containing the FORMAT statement used^ 62 to 
write the DATA records. It is initialized^- 03 "’! 06 to: 

ip, i i • • 

3 4 5 


(2HPD,[20( 


Al)) 


21 


The second word is replaced^ 4 by l 20(A6,l if 13, the third is replacedl54 
by 1 F7 . 0,| if 64. KTL words (beginning with the fifth) are re- 
placed^20-322 by |I5 , Al ,| , using the array FUSEM, which is "equivalenced" 
to FUSE(5) for convenience. 


PL0TMY 


36 




FUSEM 

FX 

FY 

I 


II 

IK 

3L 

HIM 

IM 


An array (DIMENSION is 16) "equivalenced" to FUSE(5) for convenience. 

It is initialized-*-^^ -108 bo blanks, and KTL elements are re- 
placed 320-322 by |I5,A1,| . 

The scaling parameter for the first x grid-label. It is either calcu- 
lated 2 ®® in PISTUG or supplied 232 in P(7) by the programmer if 16 . It 
is used 313 ; 3 ® 4 to compute OFX and all values of KX. If 16, it is used 234 
to compute REK. 

The scaling parameter for the first y grid-label. It is either calcu- 
lated 2 ®^ in PISTUG or supplied 244 in P(l0) by the programmer if 52 . It 
is used 313 ; 41 ® to compute OFY and all values of KY. 

The index of: 

(a) The loop 1 ®^" 1 ®® that initializes FUSEM. 

(b) The loop 21 ®" 2 ^® that computes the total number of points 

in Variation III. 

(c) The outer loop 221-225 that initializes the Y array in Varia- 

tion II. 

(d) The loop 320 " 322 that completes FUSE. 

(e) The loop 3 ®®" 3 ®® that searches for the minimum value in the 

X array. 

The Index of the inner loop 223-223 that initializes the Y array in 
Variation II. 

A counter initialized 422 to zero and used in the search loop 424 " 43 ® that 
locates the correct plotting character for Variations II and III. 

The Index of the loop 424 " 43 ® that locates the correct plotting character 
for Variations II and III. 

The number of values of Xj_ to be processed. It is the number- of times 
the loop 3 ®®" 3 ®® is executed to find a value of XMIN. ILIM is set 344 
equal to EPTS for Variation I, otherwise 34 ® to EPTST. It is also com- 
pared to J to cause branching452 when all points have been processed. 

The index of the loop 4 ®®" 44 ® that produces all KY f s and their associated 
plotting characters for one value of the X array. 


PLOTMY 


37 



near 


j 


K 


KB 


KBBC 


KBL 


KBW 


KDX 
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The subscript of the minimum value in the X array after the search 
loop 360 " 368 for XMIN is completed. (it is initialized 352 to 1.) It is 
used 402 to locate the value or values of Yj_ associated with the X-^. For 
Variations II and III, BODS’ is used in the search loop 424 " 430 that lo- 
cates the correct plotting character. BODS’ is written out 432 in an error 
message if the plotting character is not found* 

A point counter. It is initialized 337 to 1, increased 454 h y 1, and 
tested 452 for equality to ILIM to cause branching when all X-j/s have been 
processed. 

An array name, the third argument in the call list of this subroutine. 

The first three elements are always supplied by the programmer: the 

first is KOBE (see I.-B.); the second is KN (the number of curves); the 
third is a number of points to be plotted, as specified by the Variation 
selected (ref. 1, p. 11 ) . Succeeding even-numbered elements of K are 
the plotting characters, supplied by the programmer if 1 or by this, 
subroutine 422 " 1 " 8 e4 sec l- Succeeding odd-numbered elements are: 

(a) Blank if Variation I is selected. 

(b) Set 4222 equal to K(3) if Variation II is selected. 

(c) Supplied by the programmer if Variation III is selected and 
used 248 " 248 to compute EPTST. 

The output buffer area (or buffer) in which BATA records are assembled. 

KB is "equivalenced" to B. The name KB is used to move 384 ; 4481 ; 442 KX, 

KY, and the plotting character, KPCI. 

The buffer block counter. It is initialized 330 to 1, stepped 458 by 
1 each time the point counter J is stepped, used 488 to compute KBL 
(the length of the output array to be written from the buffer), and 
reset 488 to zero each time the buffer is emptied. It is also tested 470 
against zero as a simple way of recognizing the special case when EBR 
(the number of blocks in a record) is a factor of HIM (the total number 
of points). 

The' length of the array to be written 482 out of the buffer. It is the 
product 480 of KW (the number of words in a buffer block) and KBBC (the 
current total in the buffer block counter) . 

The buffer word counter. It is used as a subscript to move words into 
the buffer. Initialized 338 to 1, it is stepped 374 j 384 > 396,412,444, 44:9 
by 1 each time a word is sent to the buffer, compared 450 to EBW to see 
if the buffer is filled, and reset 484 to 1 when the buffer is emptied. 

A fixed-point name 282 for DX. It is tested 284 against zero to identify 
error message from PISTUG. 
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KDY 

KEXIT 

KHARX 

KHARY 

KK 

KLADD 

KLAST 

KLG 

KN 


KPC 


A fixed-point name 282 for DY. It is tested 284 against zero to identify 
error message from PISTUG. 

A two-way switch variable. Initialized 448 to 1, its normal sequencing 
permits the problem to continue. When all points have been calculated, 
it is set 4 ^ 4 to branch 488 to the RETURN of PLOTMY. 

An integer representing the range of the X array. It is calculated by 
PISTUG 288 or by PLOTMY 288 if L5. It is used 882 to compute KPWRX if 64 . 

An integer representing the range of the Y array. It is calculated 288 
by PIS TUG for use by PLOTXY and Is never used by PLOTMY. 

The fixed-point name for KODE. It is moved 488 from K(l) to permit 
"equivalencing" it to BKK for option branching. As KK, it is tested 448 
for zero. 

The number of characters, exclusive of KY’s and plotting characters, In 
one block of a DATA record. It is initialized 448 to 4 (three for KX, one 
for the delimiter), increased by six 488 if 8_, and by seven 488 if 64*. 

The total is then used 888 to compute KLC (the total number of characters 
in one block of a DATA record) . 

The subscript of the second last position of the K array. It is com- 
puted 848 and used as the limit of the search loop 424 ” 488 that locates 
the address of the correct plotting character for Variations II and III. 

The total number of characters in one block of a DATA record. It is the 
sum of KLADD and four* times KTL. It is computed 888 and used 882 to com- 
pute NBR, the total number of blocks possible in one DATA record (of no 
more than 132 characters) . 

The number of curves to be plotted. It is always supplied 484 by the 
programmer in K(2), and: 

(a) It is used 848 to compute KLAST. 

(b) It replaces 282 KTL for Variation I. 

(c) It is used 228 to compute KTIMES for Variation I. 

(d) It is used 228 to compute NPTST for Variations I and II. 

(e) It is used 248 to compute NPTST for Variation III. 

An array from which the standard plotting characters *, +, 0, X, =*, 0 , 
are moved 422 e4 sec L Into the K array if not 1. It is equivalenced to A. 
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KPCI The plotting character to he moved to the buffer. It is determined 

either by a search loopf 24-430 establishes the subscript of its 

position in the K array, or by direct calculation 440 of that subscript 
(possible only for Variation I) . 

KPWRX The value of 3 - KHARX. If 64, it is computed, 302 used 304 to compute 

F10X, and -written out 300 in a special heading. 

KSW1 A two-way switch variable. Initialized 114 to 1 and set 204 to 2 for 

Variation I, it is used 230 * 342 j 420 to control: 

(a) The value of NSX when PISTUG is called 200 to scale the X array. 

(b) The value of ILIM. 

(c) The method of locating the correct plotting character, KPCI. 

KSW3 A two-way switch variable. It is Initialized 110 to 1, and set 220 t'o 2 

for Variation II. It is used to control 2 ^ 0 the value of NSY when PISTUG 
is called 200 to scale the Y array,, 

KSW8 A two-way switch variable. It is initialized 113 to 1 and set 102 to 2 if 

8_. It is used 3 ^ 0 to branch to move the left-side field characters from 
the P array to the buffer. 

KSW64 A two-way switch variable. It is initialized 112 to 1, and set 132 to 2 
if 64. It is used, to branch 300 to the preparation and output of a 
special heading for the left-side field, and also to branch 300 to move 
the value of an X^ to the buffer. 

KSX One of three scaling parameters for the x scale. It is either calcu- 
lated 200 by PISTUG or supplied 231 in P(6) by the programmer if 16 . It 
is used 310 to 'compute XMF and is written out 310 in the INFO record. It 
is the number of places the decimal point must be moved to the right in 
the values of ODX and OFX by the 1401 subroutine PLOT. 

KSY One of three scaling parameters for the y scale. It is either calcu- 
lated 200 by PISTUG or supplied 242 in P(9) by the programmer if 32_. It 
is used 31 ^ to compute YMF and is written out 310 in the INFO record. It 
is the number of places the decimal point must be moved to the right in 
the values of ODY and OFY by the 1401 subroutine PLOT. 

KTIMES The value of KN-1. It is computed 220 and used only for Variation II as 
the limit of the loop 221 ” 223 that initializes the Y array. 
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KTL 


KW 


The number of Yq r s (elements in the Y array) associated with each element 
in the X array. It is initialized 111 to 1 and replaced 232 by the value 
of KN only for Variation I. It Is the number of times the loop 433 " 443 
Is executed for each Xq . In each execution one value of KY is found and 
moved to the buffer with Its associated plotting character. It is the 
total number of executions of the loop^t>-322 that completes the FORMAT 
statement FUSE. KTL Is also used 330 ; 334 to compute KLC and KW. 

The number of buffer words required to store KLC characters. It is the 
sum of KWADD and two times KTL. It is computed 334 and used 333 ; 433 to 
compute KBW and KBL. 


KWADD The number of buffer words , exclusive of those required to store KY's 

and plotting characters, needed to store one block. It is initialized 11 ? 
to 2 (one for KX and one for the delimiter), and Increased by one 139 if 
8^ and by one 139 if 64 . The total is then used 334 to compute KW. 

L A subscript computed 221 for use In the Inner loop 223 " 223 that initializes 

the Y array for Variation II. 

LL A subscript computed 432 for use in the loop 40 * 9 " 443 that produces the KY*s 

and their associated plotting characters. 

M The frequency of the x grid-lines. It is set 113+1 to 10, and, If 2 _, it 

is replaced 1 ? 2 by the value supplied In P(3) by the programmer. It is 
written out 313 in the INFO record. 

MM A value computed 222 for use in the loop 223-223 that initializes the Y 

array for Variation II. 


RBR The maximum number of blocks that can be written in one DATA record (of 
no more than 132 characters). It is computed 332 by dividing 130 (the 
number of characters in a DATA record following the "PD" required by 
1401 subroutine PLOT) by KLC (the length of a block) . It is used 33 * to 
compute EBW (the length of the buffer). 


EBW The number of buffer words required to hold the total number of charac- 
ters to be written out in one DATA record. It is computed 333 and used 430 
to test KBW (the buffer word counter) to tell if the buffer Is full. 


EF1 This is the number of characters within a block preceding the three 

character^ of KX. It is initialized 119 to zero, and is increased 133 by 
seven if £, and by six 133 if 64. The total is written out 313 In the 
UFO record. 
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NPTS 


NPTST 


NSX 


NSY 


NY 


ODX 


ODY 


OFX 


OFY 


P 


The number of values in the X array for Variation I or in the Y array for 
Variation II. It is supplied"-® 2 £,y the programmer in K(2). For Varia- 
tion I, it establishes 3 ^ the value of ILIM and is the length 232 - Q f the 
X array sent to PISTUG 2 ®® for scaling. For Variation II, it is 
stored"- 222 in K(5), K(7), . . . K(2*KN + l) , used 222 * 223 to initialize 
the Y array, used 22 ® to compute NPTST, and is the length 272 D f the Y 
array sent to PISTUG 2 ®® for scaling. 

The number of points in the X array for Variations II and III. It is 
computed 22 ® by multiplying KN (the number of curves) by NPTS for Varia- 
tion II, and by summing 2 14-218 the numbers of points for each curve 
supplied in K(3), K(5), . . . K(2*KN + l) by the programmer for Varia- 
tion III. For Variation II, it is the length 27 ® of the Y array sent to 
PISTUG for scaling; for both II and III, it is the length 23 ® of the 
X array sent to PISTUG for scaling, as well as the value 3 ^® of HIM. 

The length of the array sent 2 ®® to PISTUG to have the three scaling pa- 
rameters for the x scale calculated. For Variation I it is equal 232 to 

NPTS and for Variations II and III it is equal 23 ® to NPTST. 

The length of the array sent 2 ®® to PISTUG to have the three scaling pa- 

rameters for the y scale calculated. For Variation II it is equal 272 to 
NPTS and for Variations I and III it is equal 27 ® to NPTST. 

The frequency of the y grid-lines. It is set"-"-® +2 to 10, and if 4, it 
is replaced"-® 2 by the value supplied in P(4) by the programmer. It is 
written out 3 "-® in the INFO record. 

A number written out 3 "-® representing the x scale-factor. .It is corn- 
put ed 3 ^ 2 by multiplying the scaling parameter EX, (always an integer) by 
10"® to put it in the form required for the INFO record. 

A number written out 3 -'-® representing the y scale-factor. It is com- 
puted 3 "-^ by multiplying the scaling parameter DY, (always an integer) by 
10"® to put it in the form required for the INFO record. 

A number written out 3 "-® representing the x-scale starting -value. It is 
computed 3 "- 3 by multiplying the scaling parameter FX, (always an integer) 
by 10"® to put it in the form required for the INFO record. 

A number written out 3 "-® representing the y-scale starting-value. It is 
computed 3 "- 3 by multiplying the scaling parameter FY, (always an integer) 
by 10”® to put it in the form required for the INFO record. 

An array name, the fourth argument in the call list of this subroutine. 

The first element specifies the Variation selected; other elements are 
supplied by the programmer as required by the options chosen. 
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PFX An intermediate value found 233 and used 23 ? j_ n the calculation of KHARX. 

RD A constant established 313 equal to 5.E-7 and used 312 " 313 for rounding. 

RDX The true value of the x scale-factor. If 16 , it is computed 233 and 
used 236 to compute KHARX. 

REX The true value of the x scale starting-value. If 1(3, It Is computed 234 
and used 233 ^ 23 ? to compute KHARX. 

SFT A constant established 311 to equal l.E-6 and used 312-313 f or altering 
the scaling parameters to be written out 313 in the IHFO record. 

X An array name, the first argument In the call list of this subroutine. 

This array contains the values of the variable to be plotted on the 
vertical scale. It is sent 233 to PISTUG for scaling unless 16\ Its 
current minimum value is found (and held as XMIN) each time the search 
loop 333 " 333 is executed. That element of the array is then replaced 333 
by a large number (lO 13 ) preceding the search for a new minimum. 

XMAX An intermediate value found 23 ? during the calculation of KHARX. 

XME An intermediate value in the calculation of KX. It is computed 313 once 

and used 333 to compute each KX. 

XMIH Initialized 333 to a large number, this holds 334 the current minimum 

value of the X array at the end of the search loop. 333 ” 333 It is used 333 
to compute KX, If 64, it is multiplied 332 by F10X before being moved to 
the buffer to prepare it for printing in the left-side field. 

Y An array name, the second argument in the call list of this subroutine. 

It contains the values of the variable to be plotted on the horizontal 
scale. It is sent 233 to PISIUG for scaling unless 52 . When Variation II 
is selected, the Y array is initialized 221-223 by filling it with repe- 
titions of the given set of values. 

YMF An intermediate value in the calculation of KY. It is computed 31 ? once 
and used 433 to compute each KY. 

YS The value from the Y array currently being transformed 434 " 433 into a KY. 
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PISTUG - (FORTRAN II) 

This routine is internal to the OPP system and is used by PLOTXY and PLOTMY. 
The calling statement is: 

CALL PISTUG (KXY, NPTS, ARRAY, DELTA, FIRST, KSHIFT, KHAR) 

PISTUG accepts as input; an array in any order, its length (NPTS), and a 
switch setting (KXY, which indicates whether scaling is to be done for the x or 
y scale and whether or not the array is in order). It returns as arguments: 
suitably chosen scaling parameters (FIRST, DELTA, and KSHIFT, which together 
define a starting -value and a scale-factor) , and an integer (KHAR, which is re- 
quired in the calling routine to create a special heading when option 64 is 
used) . 

PISTUG uses the minimum and maximum values of the array to compute the total 
range. When scaling parameters are being computed for the y scale, 101 print 
positions are available to cover the range. For the x scale, an arbitrarily 
(see appendix B, V.) chosen number of lines is assigned. In either case, the 
scale factor computed is always a value of DXlO n . D is a member of the set 
[2, 2.5, 5, 10]. 

If the array is such that the field width allotted for the printed grid- 
labels (9 characters) will be exceeded (see appendix B, VI.), PISTUG writes an 
error message, sets a signal to terminate the plot, and returns control to the 
calling routine. 
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PISTUG 


SUBROUTINE PISTUG(KXY,NPTS,ARRAY,DELTA,FIRST,KSHIFT,KHAR) 


01 


SUBROUTINE PISTUG( KXY,NPTS, ARi 
C 

DIMENSION ARRAY ( 1 ) 

C 

100 N=NPTS 
102 K=KXY 
C 

104 X1=ARRAY ( 1 ) 

106 XN=ARRAY ( N ) 

107 GO T0( 302 >302, 110) , K 
C 

302 00 306 J = 1,N 
304 XI = MIN1F( X1,ARRAY( J ) ) 

306 XN = MAXlFtXN, ARRAY! J)) 

C 

110 GO TO (112,116, 112), K 
C 

112 NP=N/35+l 

113 TLIN=FL0ATF(NP)«55. 

114 GO TO 118 
C 

116 TLIN=1G1. 

C 

118 C3 = XN-Xl 
C 

120 C4=C3/TL IN 

122 KCOX = X INTF { LOGF { ABSF ( C4 ) ) / 2 
C 
C 

124 IF(25-KCUX*KCDX) 900,140,140 
C 

140 TENCDX=10.**KCDX 
142 C5=C4/ 1 ENCOX 
144 AC5=ABiF { C5 ) 

C 

150 I F ( 2. 5-AC5 ) 154,204,152 
152 IF(2.0-AC5) 204,200,200 
154 I F { 5. Q-AC5 ) 212,208,208 
C 

200 D=2. 

202 GO TO 214 
204 0=2.5 
206 GO TO 214 
208 0=5. 

210 GO TO 214 
212 0 = 10 . 

C 

214 DX=SIGNF(D,C5)*TENC0X 
C 
C 

216 GO TO (220, 230, 220), K 
C 

220 IF( (2.*C3/DX-TLIN)/TLlN-.2)222, 
C 

222 TL IN=TL IN+35. 

224 GO TO 120 


, DELTA, FIRST, K SHIFT , KHAR) 
RENAME INPUT 

FIND MIN AND MAX OF ARRAY 


TL IN = F(NPTS) OR 101 


COMPUTE TOTAL RANGE 
CALCULATE MEASURE UF SUE AND TEST 
302585 +40.) -40 

TO ERROR EXIT IF N.G. 

CALC MAX POSSIBLE DELTA 

CHOOSE BEST VALUE OF 0 


CALC SCALE FACTOR 

IF SCALING X , PLOT ENDS TOO NEAR TOP OF 
PAGE, CHANGE TOTAL LINES AND GO BACK 

1,230 
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SUBROUTINE P I STUG ( KXY , NP T S , ARRA Y , DE LT A , F I RS T , KSH I F T , KHAR ) 

C CALC STARTING VALUE 

230 AQX=AB;>F(DX) 

232 INT = AB^F(xn/ADX 
234 6EL0W=ADX*FLQATF< INT) 

236 IF(DX*Xl+l.E-20) 244,240,240 

240 F=6EL0w 

242 GO TO 246 

244 F=BELQW+ADX 

246 F X = S 1 GNF ( F , X 1 ) 

C IF SCALE VALUES INTEGRAL, GO TO 270 

250 IF(KCOX) 256,256,252 
252 KSX=6 
254 GO TO 270 

C IF NOT , ALTER AND CALC SHIFT CODE 

256 KSFT=XABsF(KCDX) +1 
258 TEMP=1U.**KSFT 
260 DX^DX* T EMP 
262 FX=FX*ItMP 
264 KSX = 6-KSFT 

C RENAME OUTPUT AND CALC KHAR 

270 DELTA = INTF 1 DX+S IGNF { . 5 , DX ) ) 

272 FIRST = INTF(FX+SIGNF( .5,FX) ) 

274 KSH I F T=KSX 

276 KHAR = X I NTF ( LOGF ( MA X IF { AB SF ( XI ) , ABSF { XN )) ) /2 . 30258 5 +40,)-40 
280 RETURN 

C 

C 

C ERROR MESSAGE 

900 WRITE OUTPUT TAPE 6,500,K, KCDX 

500 FORMAT { 30HPL PLOT HALTED BY ARRAY NUMBER , I 2 , 26H VALUES OUT OF RA 
500INGE. KCO= ,14) 

902 DX=0 • 

904 GO TO 270 

END ( I, 0,0, 0,0, 0,0, 0,0, I ,0,0,0, 0,0) 
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SUBROUTINE P I STUG ( KXY , NPTS , ARRAY, DELT A , F I RSI , KSH I FT , KHAR ) 


STORAGE NOT USED BY PROGRAM 


DEC OCT DEC OCT 

375 00567 32561 77461 


STORAGE LOCATIONS FOR VARIABLES NOT APPEARING IN COMMON, DIMENSION, OR EQUIVALENCE STATEMENT 



DEC 

UCT 


DEC 

OCT 


DEC 

OCT 


DEC 

OCT 


DEC 

OCT 

ACS 

374 

00566 

ADX 

373 

00565 

BELOW 

372 

00564 

C3 

371 

00563 

C4 

370 

00562 

C5 

369 

00561 

D 

368 

00560 

DX 

367 

00557 

F 

366 

00556 

FX 

365 

00555 

INT 

364 

00554 

KCDX 

363 

00553 

K 

362 

00552 

KSFT 

361 

00551 

KSX 

360 

00550 

NP 

359 

00547 

N 

358 

00546 

TEMP 

357 

00545 

TENCDX 

356 

00544 

TLIN 

355 

00543 

XI 

354 

00542 

XN 

353 

00541 













SYMBOLS AND LOCATIONS FOR SOURCE 

PROGRAM 

FORMAT STATEMENTS 





EFN 

LOC 


EFN 

LOC 


EFN 

LOC 


EFN 

LOC 


EFN 

LOC 

8 ) FK 

500 

00531 
















LOCATIONS 

FOR OTHER 

SYMBOLS 

NOT APPEARING 

IN SOURCE 

PRUGRAM 





DEC 

OCT 


DEC 

OCT 


DEC 

UCT 


DEC 

OCT 


DEC 

OCT 

1) 

348 

00534 

2) 

309 

00465 

3) 

314 

00472 

6 J 

327 

00507 

9) 

346 

00532 

CIGl 

351 

00537 

C ) G2 

352 

00540 

E 1 4 

64 

00100 

E I 6 

85 

00125 

E)B 

140 

00214 

E) 0 

150 

00226 

E) I 

178 

00262 

E ) K 

213 

00325 

E ) 0 

232 

00350 

E)Q 

296 

00450 

E ) 10D 

150 

00226 

E ) 20 J 

184 

00270 















LOCATIONS 

OF NAMES IN TRANSFER 

VECTOR 







DEC 

OCT 


DEC 

OCT 


DEC 

OCT 


DEC 

OCT 


DEC 

OCT 

EXP ( 2 

1 

00001 

LOG 

0 

00000 

(FIL) 

3 

00003 

(STH) 

2 

00002 





ENTRY POINTS TO SUBROUTINES NOT OUTPUT FROM LIBRARY 
EXP { 2 LOG { F I L ) (STH) 

EXTERNAL FORMULA NUMBERS WITH CORRESPONDING INTERNAL FORMULA NUMBERS AND OCTAL LOCATIONS 


EFN 

IFN 

LOC 

EFN 

IFN 

LOC 

EFN 

IFN 

LUC 

EFN 

IFN 

LOC 

EFN 

IFN 

LOC 

100 

5 

00042 

102 

6 

00045 

104 

7 

00050 

106 

8 

00052 

107 

9 

00054 

302 

10 

00057 

304 

11 

00062 

306 

12 

00070 

110 

13 

00102 

112 

14 

00105 

113 

15 

00115 

114 

16 

00124 

116 

17 

00127 

118 

18 

00131 

120 

19 

00134 

122 

20 

00137 

124 

21 

00156 

140 

22 

00166 

142 

23 

00174 

144 

24 

00177 

150 

25 

00202 

152 

26 

00207 

154 

27 

00216 

200 

28 

00223 

202 

29 

00225 

204 

30 

00230 

206 

31 

00232 

208 

32 

00233 

210 

33 

00235 

212 

34 

00236 

214 

35 

00240 

216 

36 

00246 

220 

37 

00251 

222 

38 

00264 

224 

39 

00267 

230 

40 

00272 

232 

41 

00275 

234 

42 

00307 

236 

43 

00317 

240 

44 

00327 

242 

45 

00331 

244 

46 

00332 

246 

47 

00335 

250 

48 

00341 

252 

49 

00345 

254 

50 

00347 

256 

51 

00352 

258 

52 

00356 

260 

53 

00364 

262 

54 

00367 

264 

55 

00372 

270 

56 

00375 

272 

57 

00404 

274 

58 

00413 

276 

59 

00415 

280 

60 

00444 

900 

62 

00452 

902 

64 

00462 

904 

65 

00464 





-r 



100-106 


PISTUG 


1 ^ 

DX = 0 


1 *■ 

(Error signal) 




Compete KSX 
and complete 
DX,FX 


r V 

r 246 1 

Get 

sign for 


FX 

250 1 

\ 

KCDX / 

7 

A 

P 270 
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PIS TOG - DICTIONARY 


Superscripts are statement numbers in the FORTRAN listing (p. 45) . 

The notation £ means option z is being used. 

AC5 When the maximum possible scale-factor , C4, is expressed as x.xxxx multi- 
plied by 10 n , AC5 is 144 the absolute value of the x.xxxx factor. It is 
used 150 " 154 to select D, the next larger value from the set of permis- 
sible values [2, 2.5, 5, 10]. 

ADX The absolute value of DX. It is used 252 " 248 in the calculation of the 
scaling parameter FX. 

ARRAY Input to PISTOG, the third argument in the call list. Its minimum and 
maximum values are found in a loop 502 ~ 508 if they are not already known 

(as is the case when PISTOG is scaling an X array for PLOTXY) . 

BELOW An intermediate step in the calculation 230 " 248 Q f px. 

C3 The total range of the array being scaled. It is computed 118 from the 

minimum and maximum values of the array and used 120 to compute C4 (the 
maximum possible scale-factor) . When an X array is being scaled, C3 is 
used 220 to decide if TLIN should be modified. 

C4 The value of the maximum possible scale-factor. It is found 120 by 

dividing the range, C3, by TLIN and used 122 to compute KCDX. 

C5 When the maximum possible scale-factor C4 is expressed as x.xxxx multi- 

plied by 10 n , C5 is the x.xxxx factor. It is computed 142 and used 144 to 
find the value of ACS. 

D This is found 150 " 212 as the smallest member of the set [2, 2.5, 5, 10] 

that is larger than ACS. It establishes 214 the first value of the 
scaling parameter DX. 

DELTA Output from PISTOG, the fourth argument in the call list. It is the 

scaling parameter that represents the scale-factor. It is established 
as EX and renamed 270 to avoid excessive address modification. 

DX The scaling parameter that represents the scale-factor. It is origi- 

nally 214 the true value of the scale-factor; if KCDX < 0, DX is multi- 
plied 280 by a suitably chosen power of 10 and is stilT named EX. DX is 
set 902 to zero as an error signal to the calling routine. 

F An intermediate step in the calculation 230-246 of FX. 
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FIRST 


FX 

HOT 

K 


KCDX 

KHAR 

KSFT 

KSHIFT 


Output from PISTUG, the fifth argument in the call list. It Is the 
scaling parameter that represents the star ting -value. It is established 
as FX and renamed 212 to avoid excessive address modification. 

The scaling parameter that represents the starting -value. It is origi- 
nal 1y 24 ® the true value with which the first grid-line will be labeled; 

If KCDX <0, FX is multiplied 2 ® 2 by a suitably chosen power of 10 and is 
still named FX. 

An intermediate step in the calculation 252 of FX. 

A switch variable, supplied 102 (as KXY) by the calling routine. It will 
have the value: 

(a) 1 when CALL is to scale X from PLOTMY. 

(b) 2 when CALL is to scale Y from either PLOTXY or PLOTMY. 

(c) 3 when CALL is to scale X from PLOTXY. 

It is used 101 * 110 * 21 ® to control branching: 

(a) To find the minimum and maximum of the array. 

(b) To determine TLIN, 

(c) To reevaluate TLIN if necessary. 

This is the characteristic of log-^o C4. It is computed 122 and tested 124 
to see if it is equal to or between -5 and +5, and if it Is not, an error 
message is written 900 and an error signal (DX=0) returned to the calling 
routine. It is also used 140 to compute TENCDX. If KCDX is zero or nega- 
tive 2 ® 0 , the scaling parameters FX and DX are multiplied by 10 K°IOT. 

(KSFT = | KCDX | + 1) 

Output from PISTUG, the seventh argument of the call list. It is calcu- 
lated 21 ® by finding the characteristic of the log-^Q of either the maxi- 
or the minimum of the array, whichever has the greater absolute value. 

An intermediate value calculated 250 whenever KCDX < 0. It is used 25 ® as 
the power of 10 by which the scaling parameters FX and DX are multiplied. 

Output from PISTUG, the sixth argument In the call list. This scaling 
parameter is the decimal-point-shift code applicable to the other two 
scaling parameters DELTA and FIRST, which it accompanies. It is com- 
puted as KSX and renamed 214 to avoid excessive address modification. 


PISTUG 


50 


KSX 


The scaling -parameter representing the decimal -point -shift code for DX 
and EX. If KCDX is greater than zero, KSX is set 232 equal to six* If 
KCDX is less than or equal to zero, KSX is set 238 " 284 equal to 
5 - |kcdx|. 

KXY Input to PISTUG, the first argument in the call list* It is renamed 1 * 32 K 
to avoid excessive address modification. 

N The number of points in the array being scaled. It is supplied 18 ** 3 by 

the calling routine (as HPTS) in the call list. It is used as the limit 
of the loop 3 * 32 " 3 * 38 that locates the minimum and maximum values of the 
array and in the original calculation 112 " 114 of TLIK when scaling an 
X array. 

HP An intermediate value computed 112 in the calculation of TLITL 

HPTS Input to PISTUG, second argument in the call list. It is renamed 18 ** 3 to 
avoid excessive address modification. 

TEMP Temporary storage. 

TEKCDX The value of ioPCDX^ it is computed 148 ^ and used to compute C5 142 and 
DX. 214 

TLIK When PISTUG is scaling a Y array, TLIK is set 118 equal to 101. When an 
X array is being scaled, TLIU is found 113 first as an empirical function 
of the length of the array, later it may be modified. 220-222 (See ap- 
pendix B, V.) In either case, TLHF is used 120 " 122 to compute KCDX. 

XI The minimum or maximum value of the array being scaled. It is set 18 * 4 

equal to the first element of the array when PISTUG is scaling an X array 
for PLOTXY, otherwise it contains 38 * 4 the minimum value of the array 
after the search loop 3 * 32 " 38 * 8 is completed. It Is used 118 to compute C3, 
the range of the array, and enters 232 j 238 .> 248 several times into the 
calculation of EX, the scaling parameter representing the star ting -value. 
It is also used 218 to compute KHAR. 

XK The minimum or maximum value of the array being scaled. It is set 18 * 8 

equal to the last element of the array when PISTUG is scaling an X array 
for PLOTXY, otherwise it contains 3 * 38 the maximum value of the array 
after the search loop 38 * 2 * -38 * 8 ± s completed. It is used 118 to compute the 
range of the array and in the calculation 2 ^ 8 of KHAR. 
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SORTXY - (FORTRAN II) 


The calling statement is: 

CALL SORTXY (V, W, NPTS) 

This subroutine rearranges the NPTS elements of the V array in order of in- 
creasing size. The elements of the W array are moved to maintain the original 
pair -relation; that is, if the fifth element of the V array is moved to the 
first position of V, the fifth element of W is moved to the first position of W. 
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SORTXY 


SUBROUTINE SORTXYI X, Y, NPTS ) 


SUBROUTINE SORTXY( X, Y,NPTS ) 
DIMENSION X(IQO) »Y( IOO) 

IOO N=NPTS 
102 NN=N-i 
104 DO 140 KT=1,NN 
XMIN=XIKT) 

JAD*KT 

JKL=KT+1 

112 DO 120 JK= JKL ,N 

114 IF ( XM IN-X ( JK ) ) 120,120,116 

116 XMIN^XiJK) 

118 JAD-JK 
120 CONTINUE 
l 22 YMIN=Y( JAD) 

XUAD)* X ( KT ) 

Y ( JAD ) = Y ( KT) 

X ( KT) * XMIN 
Y ( KT ) ~ YMIN 
140 CONTINUE 
RETURN 

END( 1,0, 0,0, 0,0, 0,0, 0,1,0, 0,0, 0,0) 


SUBROUTINE SQRTXYI X, Y, NPTS ) PAGE 2 


STORAGE NOT USED BY PROGRAM 

DEC OCT DEC OCT 

9? 00141 32561 77461 


STORAGE LOCATIONS FOR VARIABLES NOT APPEARING IN COMMON, DIMENSION, OR EQUIVALENCE STATEMENT 



DEC 

OCT 

DEC 

OCT 


DEC OCT 


DEC OCT 


DEC 

OCT 

JAD 

96 

00140 

JKL 95 

00137 

JK 

94 00136 

KT 

93 00135 

NN 

92 

00134 

N 

91 

00133 

XMIN 90 

00132 

YMIN 

89 00131 









LOCATIONS 

FOR OTHER 

SYMBOLS 

NOT APPEARING 

IN SOURCE 

PROGRAM 





DEC 

OCT 

DEC 

OCT 


DEC OCT 


DEC OCT 


DEC 

OCT 

2) 

80 

00120 

6) 81 

00121 

9) 

87 00127 

C )G2 

88 00130 






EXTERNAL 

FORMULA NUMBERS WITH CORRESPONDING INTERNAL FORMULA NUMBERS AND OCTAL 

LOCATIONS 



EFN 

IFN 

LOC 

EFN IFN 

LOC 

EFN 

IFN LOC 

EFN 

IFN LOC 

EFN 

IFN 

LOC 

100 

4 

00031 

102 5 

00033 

104 

6 00036 

112 

10 00054 

114 

11 

00061 

116 

12 

00066 

118 13 

00070 

120 

14 00074 

122 

15 00077 

140 

20 

00111 



SCALE - (FORTRAN II) 


The calling statement is: 

CALL SCALE (EPTS, A, KRSTR) 

This subroutine finds the largest absolute value of the NPTS elements of A, 
then computes the characteristic of its log^g* If the characteristic K is 
-2 < K < 4, KRSTR is set to zero and control returns to the calling program. If 
K >4 or < -2, each element of A is multiplied by a power (KRSTR) of 10 to trans- 
form the array to suit PLOTXY and PL01MY. 

KRSTR is output to permit the user to ReSToRe the array or record how it 
has been altered. 
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SCALE 


SUBROUTINE SCALE ( NPTS, X , KRSTR ) 


SUBROUTINE SCALE t NPTS , X, KRSTR ) 

DIMENSION X(l) 

100 N=NP T S 
102 XMAX = XU) 

104 00 106 J=1,N 

106 XMAX=MAX1F{ ABSF(XMAX) ,ABSF( X( J 1 ) ) 

108 KHAR = XlNTFlL0GF(XMAX>/2. 302585+40. >-40 
116 I F ( ( 4- KHAR ) * ( KHAR+2 ) ) 120,118, 118 

118 KP WR=0 

119 GO TO 130 

120 KPWR=3-KHAR 
122 FACT=10.**KPWR 
126 00 120 J=1 » N 
128 XI J)=X(J)*FACT 
130 KRSTR=KPWR 

140 RETURN 

END (1,0, 0,0, 0,0,0, 0,0,1, 0,0, 0,0,0) 


SUBROUTINE SCALE ( NPTS , X ,KRSTR ) 


STORAGE NOT USED BY PROGRAM 

DEC OCT DEC OCT 

118 00166 32561 77461 


STORAGE LOCATIONS FOR VARIABLES NOT APPEARING IN COMMON, DIMENSION, OR EQUIVALENCE STATEMENT 


FACT 

DEC 

117 

OCT 

00165 

KHAR 

DEC 

116 

OCT 

00164 

KPWR 

DEC 

115 

OCT 

00163 

N 

DEC 

114 

OCT 

00162 

XMAX 

DEC 

113 

OCT 

00161 




LOCATIONS 

FOR OTHER 

SYMBOLS 

NOT APPEARING 

IN SOURCE 

PROGRAM 




1) 

DEC 

no 

OCT 

00156 

2) 

DEC 

93 

OCT 

00135 

3) 

DEC 

99 

OCT 

00143 

6) 

DEC 

102 

OCT 

00146 

9) 

OEC 

108 

OCT 

00154 






LOCATIONS 

OF NAMES IN 

TRANSFER 

VECTOR 






EXP ( 2 

DEC 

1 

OCT 

00001 

LOG 

DEC 

0 

OCT 

00000 


DEC 

OCT 


DEC 

OCT 


OEC 

OCT 


ENTRY POINTS TO SUBROUTINES NOT OUTPUT FROM LIBRARY 

EXP ( 2 LOG 


EXTERNAL FORMULA NUMBERS WITH CORRESPONDING INTERNAL FORMULA NUMBERS AND OCTAL LOCATIONS 


EFN 

1FN 

LOC 

EFN 

IFN 

LOC 

EFN 

IFN 

LOC 

EFN 

IFN 

LOC 

EFN 

IFN 

LOC 

100 

4 

00026 

102 

5 

00030 

104 

6 

00032 

106 

7 

00035 

108 

8 

00052 

116 

9 

00070 

118 

10 

00103 

119 

11 

00105 

120 

12 

00106 

122 

13 

00111 

126 

14 

00117 

128 

15 

00122 

130 

16 

00127 

140 

17 

00131 






PLOT - (1401 SPS) 


Control is transferred to PLOT from the Lewis 1401 tape-to-printer program 
whenever a record has been read whose first character is a P. PLOT recognizes 
as input, and processes differently, four types of records; those whose first 
two characters are: 

PT - hereinafter called TITLE records. 


PI 

n 

it 

INFO 

If 

PD 

it 

ii 

DATA 

tl 

PL 

ir 

n 

LEGEND 

If 


The result of this processing is a printed plot consisting of x and y grid- 
lines, x and y grid-labels, plotted points, and possibly printed information at 
the top, bottom, and to the left of the grid. Sample tape records and the corre- 
sponding plotted output are in appendix A. 

X. The input to this routine is any and all BCD records whose first charac- 
ter is P. The routine recognizes and processes any such record whose 
second character is T, I, D, or L. If any other character is in the second 
position of the record, the record is printed and control is returned to 
the tape-to-printer program. 

A. TITLE records. 

The remaining contents of a TITLE record (starting with the third 
character) are printed above the plot, starting in print position 1. 

The PT is not printed. The total length of a TITLE record is limited 
to 132 characters . 

B. IKFO records. 

The remaining contents of an INFO record (starting with the third 
character) are information needed by PLOT to determine the construction 
of the grid, the construction of the grid-labels, the plotting character 

or characters, and the format of the DATA records that follow. Specifi- 
cally, the PI record consists of the characters PI followed by 10 fields 

whose contents are: 

(Field width is shown in parentheses.) 

a. (l) This is the plotting character for all points unless it is a 

In this case, the plotting characters will be supplied in the 
following DATA records. 

b. (l) This digit is the decimal -point-shift code, EXSH, for XSTART 
and DX. (items d. and f. in this list.) 
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c. (l) This digit is the decimal-point-shift code, PYSH, for YSTART 
and DY. (items e. and g. in this list.) 

d. (9) The number (XSTART) in this field is in the form tO.xxxxxx and 
represents the first x grid-label. 

e. (9) The number (YSTART) in this field is in the form tO.xxxxxx and 
represents the first y grid-label. 

f. (9) The number (DX) in this field is in the form tO.xxxxxx and 
represents the x scale-factor (value of one line-space) . 

g. (9) The number (DY) in this field is in the form tO.xxxxxx and 
represents the y scale-factor (value of one print position) . 

h. (3) This field contains an integer, XGD, specifying the frequency 
of the x grid-lines. They will be printed every XGD line-spaces below 
the first. If XGD = 0, only the first and final x grid-lines will be 
printed. 

i. (3) This field contains an integer, YGD, specifying the frequency 
of the y grid-lines. They will be printed every YGD positions to the 
right of the first. If YGD = 0, only the first will be printed. 

j. (3) This is an integer, FM1, specifying the length of the first 
field in each block of information in the following DATA record. 

C. DATA records. 

The remaining contents of a DATA record starting with the third 
character are successive blocks of information needed by PLOT to position 
the plotting characters and to print the left-side field. Only complete 
blocks are permitted within a record. Within blocks, the location in 
which a plotting character is to be printed is specified by two positive 
integers, X, the number of line-spaces below the first x grid-line, and Y, 
the number of print positions to the right of the first y grid-line. 
Specifically, the format of a block is: 

(a) The first FM1 characters are the left-side field. 

(b) The next three characters are X. 

(c) (l) If the first character in the INFO record (after PI) is not a 
minus sign, the next three characters are Y and the block is complete. 

(2) If it is a minus sign, the next three characters are Y and the 
following character is the plotting character for that point j since 
additional points may fall on the same line, additional Y f s and their 
associated plotting characters may follow the first. The end of a 
block is indicated by a one-character delimiter, the letter E. 
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D. LEGETTO records. 


The remaining contents of a LEGEND record starting with the third 
character are printed below the plot starting in print position 1. The 
PL is not printed. The total length of a LEGEND record is limited to 
132 characters . 

II. The output from this routine consists of: 

a. The contents of the first TITLE record printed at the top of a new 
page. The contents of succeeding TITLE records are printed single space, 
followed by a single space. 

b. A line of y grid-labels. They are constructed using YSTAPT, DY, and 
FYSH. (See I.-B.-c., e.,g.) The first label is centered over the first 
y grid-line, and the others are centered successively around every tenth 
print position to the right. The maximum length of a label is nine print 
positions; seven digits, decimal point, and sign. 

c. A set of y grid-lines. The first is in print position 30 and the 
others are in every YGD th print position to the right (i.-B.-i.). The 

y grid-line character is "l". If a plotted point falls on the grid-line, 
the "l” is replaced by the appropriate plotting character. 

d. A set of x grid-labels. They are constructed using XSTART, DX, and 
PXSH (I. -B.-b.,d.,f. ) . The first is printed with the first x grid-line 
to the immediate left of the first y grid-line. The others are printed 
every 10 line-spaces below this until the plot is terminated. The maxi- 
mum label length is nine print positions; seven digits, decimal point, 
and sign. 

e. A set of x grid-lines. The first is printed immediately below the 

line of y grid-labels (described in II. -B.) and the others every XGD 
line-spaces below the first until the last point has been printed 
(i.-B.-h.). If it did not lie in a grid line, the plot is continued un- 
til the next x grid-label is printed and a final x grid-line is then 
printed with the label. The x grid-line character is If a plotted 

point falls on the grid-line, the Tf ~ M is replaced by the appropriate 
plotting character. 

f. Points. Plotting characters are positioned on the page using X and Y 

(I.-C.). 

g. Up to 19 characters to the left of the first y grid-line on any line 
on which a point is plotted. This is called the left-side field. 

h. A line of y grid-labels (similar to II. -b.) is printed below the 
final x grid-line. 
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i. A double space follows the last line of grid-labels before the con- 
tents of the first LEGEND record are printed. Successive LEGEND records 
are printed single space. 

III. PLOT is an SPS subroutine called by the Lewis 1401 tape-to-printer 
program that processes 7090 output tapes. An input record is expected to 
be in core, starting at RECORD, each time PLOT is called. 

Each plot is started at the top of a new page. One or more TITLE 
records are printed, a single space occurs, and a line of y grid-labels 
(centered around every 10 print positions) is printed. 

The first line of the plot is then constructed in the print area. It 
consists of the first x grid-label, the first x grid-line, the y grid-line 
character (every 10 print positions), and the plotting characters for any 
points with a line count of zero. Once all of the points whose line count 
is zero have been processed, the line is printed. Succeeding lines of the 
plot are constructed in the same manner with the following exceptions: 

1. The x grid-labels are printed only on every tenth line. 

2. The x grid-line frequency of printing is controlled by the input 
value XGD. 

3. The spacing of bhe y grid-line characters is controlled by the y grid- 
line frequency, YGD. 

When the first LEGEND record is received (indicating no more data 
points), the plot is completed in the following manner. If the line ready 
in the print area contains an x grid-line, it is printed and followed by a 
line of y grid-labels. If there is no x grid-line in the print area, then, 
after the current line is printed, the y grid-lines are continued until the 
line for the next x grid-label is reached. A final x grid-line accompanies 
this label and is followed by a line of y grid-labels. A double space 
occurs before the waiting LEGEND record is printed, and any futher LEGEND 
records are printed single space. 

After the routine has processed any one record beginning with P, con- 
trol is returned to the calling program. 
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PLOT 


1401 SUBROUTINE PLOT 


THIS ROUTINE IS WRITTEN IN AN ASSEMBLY LANGUAGE SIMILAR TO 1401 
SPS. THE ASSEMBLER IS A 7090 FORTRAN II PROGRAM WHICH HAS BEEN 
DISTRIBUTED THROUGH SHARE (NUMBER 1080). 

THIS LISTING WAS PREPARED BY A SPECIAL 1401 PROGRAM WHICH COMBINED 
COMMENTS CARDS WITH THE ASSEMBLER OUTPUT. 

NO SYMBOL TABLE DICTIONARY IS INCLUDED FOR THIS SUBROUTINE* SINCE 
THE DATA DEFINITIONS AND THE ALPHABETIC CROSS-REFERENCED LABEL 
TABLE WHICH FOLLOW SERVE THE SAME PURPOSE. 


NO. 

CC 

LABEL 

OP 

A-ADD. 

B-ADD. 

D 

------- COMMENTS ------- 

LOCATION 

1 

00 


ORG 

1000 



ALL EXITS FROM PLOT 


hd 

2 2 

07 

PEXIT 

MCW 

EDIT 

EDX 


WHEN A PLOT IS FINISHED (OR STOPPED WITH AN ERROR MESSAGE), THE 

1000 

fa 3 

01 


MCW 




EDIT WORDS MUST BE CLEARED, SWITCH FSTREC SET FOR FIRST-TITLE 

1007 

4 

07 


SW 

FSTREC 

LSTLAB 


AND SWITCH LSTLAB SET FOR NOT-FIRST LEGEND - - BEFORE CONTROL 

1008 

5 

04 

START 

B 

• 



IS RETURNED TO TTP THROUGH RETURN STORED AT *START*. 

1015 








ENTRY TO 









SUBROUTINE PLOT 


6 

04 

ENTRY 

SBR 

START + 3 



STORE RETURN AT * START*. 

1019 

7 

04 


SW 

PPGUP 



SET SKIP-TO-NEW-PAGE-SWITCH (PPGUP) FOR USE BY TTP. 

1023 

8 

08 


B 

DATA 

REC +1 

D 

IF DATA RECORD, GO TO *DATA*. IF INFO RECORD, GO TO *FORMAT*. IF 

1027 

9 

08 


B 

FORMAT 

REC +1 

I 

LEGEND RECORD, GO TO *LABEL*. IF NONE OF THESE, TREAT AS A 

1035 

10 

08 


B 

LABEL 

REC +1 

L 

TITLE RECORD. 

1043 


PLOT 


NO. CC 

LABEL OP 

A-ADD. 

B-ADD. 

0 

— 

— 

— 

— 

— 

- COMMENTS - 


— 

* 

LOCATION 












PROCESS 

TITLE 

RECORDS 



11 

01 


NOP 



FSTREC-F IRST-T ITLE SWITCH. 

1051 

12 

02 

FSTREC 

CC 



1 IF FIRST TITLE# SKIP TO NEW PAGE. 

1052 

13 

07 

TITLEP 

cw 

FSTPT 

FSTREC 

SET SWITCH FSTPT FOR FIRST-DATA AND SWITCH FSTREC FOR 

1054 







NOT-FIRST TITLE. 


14 

04 


cs 

336 


CLEAR PRINT AREA. 

1061 

15 

01 


cs 




1065 

16 

07 


MCM 

REC +2 

201 

MOVE TITLE RECORD! OM ITT I NG 1ST 2 CHARACTERS-PT-) INTO PRINT AREA. 

1066 

17 

07 


MCW 

EDIT -8 

336 

MOVE X GRID-LINE CHARACTER {-) INTO 336. 

1073 

18 

04 


w 

START 


PRINT THE TITLE RECORD AND EXIT TO TTP THRU *START*. 

1080 







PROCESS DATA RECORDS. 



19 

07 

DATA 

SbR 

A 

REC +1 


INITIALIZE INOEX A FOR STEPPING THROUGH DATA RECORD. 

1084 

20 

01 


NOP 




FSTPT-FIRST-DATA SWITCH 

1091 

21 

04 

FSTPT 

B 

ZX 



IF NOT-FIRST DATA, GO TO *ZX*. 

1092 

22 

04 


SW 

FSTPT 



IF FIRST-DATA, SET SWITCH FSTPT FOR NOT-FIRST DATA. 

1096 

23 

04 

INTX 

S 

LASTX 



CLEAR LINE COUNTER ( LASTX) TO ZERO. 

1100 

24 

07 

MIYRET 

C 

XLBLTS 

TEN 


IF IT IS TIME FOR AN X GRID-LABEL, (EVERY 10TH LINE), GO TO 

1104 

25 

05 


B 

XGL 


S 

SUBROUTINE *XGL* TO POSITION LABEL IN PRINT AREA. 

1111 

26 

07 

XLRET 

C 

YGDTST 

XGD 


IF IT IS NOT TIME FOR AN X GR ID-LINE, ( EVERY-XGD-L INES ) , 

1116 

27 

05 


& 

YGDR 


/ 

GO TO *YGDR*. 

1123 

28 

04 

YGRID 

SW 

PLZERO 



IF IT IS, PREPARE TO PUT X GRID-LINE CHARACTER INTO PRINT 

1128 

29 

07 


SBR 

YGDTST 

0 


AREA. CLEAR COUNTER YGDTST TO ZERO. FILL PRINT AREA (EXCEPT 

1132 

30 

07 

NOYGD 

MCW 

MINUS 

335 


LEFT-SIDE FIELD) WITH MINUS SIGNS. 

1139 

31 

04 

YGDR 

B 

VTGRID 



GO TO SUBROUTINE *VTGRID* TO POSITION Y GRID-LINE CHARACTERS IN 

1146 








PRINT AREA IN EVERY YGD PRINT POSITIONS. 


32 

04 

ZX 

S 

Y 



CLEAR WORKING AREA FOR X. 

1150 

33 

04 


MCW 

Y 



CLEAR WORKING AREA FOR Y. 

1154 

34 

07 

XADD 

A 

0 

A X 


MOVE X IN FROM DATA RECORD. 

1158 

35 

07 

YADD 

A 

0 

A Y 


MOVE Y IN FROM DATA RECORD. 

1165 

36 

07 

XSTOl 

C 

X 

LASTX 


IF X=LASTX( LINE COUNTER), POINT IS TO BE PRINTED ON THIS LINE, SO 

1172 

37 

05 

TNXTX 

B 

PLOT 


s 

GO TO PLOT. 

1179 

38 

07 


S 

LASTX 

X 


IF NOT, FORM (X-LASTX). IF THIS VALUE IS NEGATIVE, GO TO 

1184 

39 

08 

DXTST 

BWZ 

STOP 

X 

K 

ERROR MESSAGE VIA *STOP*. 

1191 

40 

07 


A 

LASTX 

X 


IF OK, RESTORE VALUE OF X. 

1199 

41 

01 


W 




PRINT A LINE. 

1206 

42 

04 


CS 

332 



CLEAR PRINT AREA AND INITIALISE FOR NEXT LINE. 

1207 

43 

01 


cs 





1211 

44 

07 


A 

ONE 

XLBLTS 


STEP COUNTER FOR X GRID-LABELS (XLBLTS) BY 1. 

1212 

45 

01 


A 




STEP LINE COUNTER (LASTX) BY 1. 

1219 

46 

01 


A 




STEP COUNTER FOR X GRID-LINES (YGDTST) BY 1. 

1220 

47 

07 


A 

DX 

XSTART 


STEP VALUE OF XSTART (USED FOR X GRID-LABELS) BY DX. 

1221 

48 

04 

CHGSIG 

B 

MIYRET 



GO TO «MI YRET*. 

1228 


PLOT 


NO. 

CC 

LABEL 

OP 

A-ADD. 

B-ADD. 


D 

------- COMMENTS ------- 

LOCATION 

49 

04 

PLOT 

SW 

1 

A 



SET WORD-MARK IN DATA RECORD TO END MOVE. 

1232 

50 

07 

PLOTA 

MCW 

0 

A 219 



MOVE LEFT-SIDE FIELD TO PRINT AREA. 

1236 

51 

04 


CW 

1 

A 



CLEAR WORD-MARK. 

1243 

52 

07 

YST01 

MN 

Y 

B 



STORE Y IN INDEX B (OMITTING SIGN). 

1247 

53 

01 


MCW 






1254 

54 

07 

YST02 

C 

B 

YLAST 



IF Y IS TOO BIG (GREATER THAN 103), GO TO *TQBIG*. 

1255 

55 

05 


B 

TOBIG 



T 


1262 

56 

01 


NOP 





MCBRA -MULT I PL E-CUR VE SWITCH. 

1267 

57 

04 

MCBRA 

B 

MANYY 




IF MULTIPLE CURVES, GO TO *MANYY*. 

1268 

58 

07 


MCW 

PLCHAR 

PLZERO 

B 


IF SINGLE CURVE, MOVE PLOT CHARACTER INTO PRINT AREA B 
POSITIONS TO THE RIGHT OF THE FIRST Y GRID-LINE. 

1272 

59 

07 

PLGTBK 

A 

TWONE 

A 



STEP A TO NEXT BLOCK OF DATA RECORD. 

1279 

60 

08 

ONWARD 

B 

START 

6 

A 

* 

END OF RECORD TEST. 

1286 

61 

01 


B 





IF ONE OF THE NEXT SIX CHARACTERS IN THE DATA RECORD 

1294 

62 

01 

• 

B 





IS A GROUP MARK, GO TO *ST ART*. IF NOT, GO 

1295 

63 

01 


B 





TO *ZX* TO PROCESS NEXT BLOCK. 

1296 

64 

01 


B 






1297 

65 

01 


B 






1298 

66 

04 


B 

ZX 





1299 

67 

01 

MANYY 

NOP 





MAN YY=F IRST-Y-IN-BLGCK SWITCH. 

1303 

68 

04 


B 

MODA 




IF THIS IS FIRST Y, BRANCH TO *MOOA* AND RETURN. 

1304 

69 

04 

MUPLC 

SW 

PLZERO 

B 



SET WORD-MARK FOR PLOT CHARACTER AND MOVE PLOT CHARACTER INTO 

1308 

70 

07 


MCW 

1 

A PLZERO 

B 


PRINT AREA B POSITIONS TO THE RIGHT OF THE FIRST Y GRID-LINE. 

1312 

71 

07 

UPA 

SBR 

A 

1 

A 


STEP INDEX A BY 1. 

1319 

72 

08 

LTYTST 

BBE 

ENOTST 

1 

A 

+ 

IF NEXT CHARACTER IS THE DELIMITER, GO TO *ENDTST*. 

1326 

73 

07 


SBR 

A 

3 

A 


IF NOT, STEP INDEX A FOR ADDRESS OF NEXT Y. 

1334 

74 

04 


S 

Y 




CLEAR WORKING STORAGE Y. 

1341 

75 

07 


A 

0 

A Y 



MOVE NEXT Y IN FROM DATA RECORD. 

1345 

76 

04 


B 

YST01 




GO BACK TO *YSTQ1 * • 

1352 

77 

07 

ENDTST 

SBR 

A 

1 

A 


STEP INDEX A TO NEXT BLOCK OF DATA RECORD. 

1356 

78 

04 


SW 

MANYY +1 




SET SWITCH MANYY FOR F1RST-Y. 

1363 

79 

04 


B 

ONWARD 




GO TO *ONWARD*. 

1367 

80 

04 

MODA 

CW 

MANYY +1 




SET SWITCH MANYY FOR NOT-FIRST-Y. 

1371 

81 

07 


A 

TWONE 

A 



STEP INDEX A TO LOCATE FIRST PLOT CHARACTER IN DATA BLOCK. 

1375 

82 

04 


B 

MUPLC 




GO TO *MUPLC*. 

1382 

83 

07 

T0B1G 

MCW 

RM 

332 



PUT A RECORD MARK AT RIGHT END OF PRINT AREA TO SHOW Y TOO BIG. 

1386 

84 

08 


BWZ 

UPA 

MCBRA 


1 

IF SWITCH MCBRA IS SET FOR MULTIPLE CURVES, GO TO *UPA*. IF FOR 

1393 

85 

04 


B 

PLQTBK 




SINGLE CURVE, GO TO *PLOTBK*. 

1401 


CD 



LOCATION 


CD 


HU 

It* 


§ 


NO. 

CC 

LABEL 

OP 

A-ADD. 

B-ADD. 


D 

86 

04 

NXREC 

SBR 

ENDNR +3 




87 

04 


CS 

332 




88 

01 


CS 





89 

07 

LABOUT 

MCW 

YSTART 

TEMPI 



90 

07 


SBR 

A 

0 



91 

07 

NXTLAB 

SBR 

B 

0 



92 

07 


C 

TEMPI -1 

ZTST -1 



93 

05 


B 

EDYST 



/ 

94 

07 

YZRTST 

MN 

TEMPI 

PLZERO 

A 


95 

08 


B 

YZROK 

PLZERO 

A 

0 

96 

07 

EDYST 

LCA 

EDY 

TEMP2 



97 

07 


MCE 

TEMPI 

TEMP2 



98 

08 

TSTB 

BBE 

NONB 

TEMP2 -6 

B 

* 

99 

07 


SBR 

B 

1 

B 


100 

04 


B 

TSTB 




101 

04 

NONB 

SW 

TEMP2 -7 




102 

07 


MZ 

TEMP2 -8 

TEMP2 -7 

B 


103 

08 

UPB 

B 

MVLAB 

B 


6 

104 

07 


SBR 

A 

1 

A 


105 

07 


SBR 

B 

1 

B 


106 

07 


A 

199 

TEN 



107 

04 


B 

UPB 




108 

07 

MVLAB 

MCW 

TEMP2 

PLZERO— 2 

A 


109 

07 

YZROK • 

A 

TEN 

A 



110 

07 


SBR 

TEN 

10 



111 

07 


A 

DY 

TEMPI 



112 

07 


C 

ONETEN 

A 



113 

05 


6 

ENDNR 



U 

114 

04 


B 

NXTLAB 




115 

04 

ENDNR 

W 

• 





116 

01 

STOP 

W 




117 

04 


CS 

332 



118 

01 


CS 




119 

07 


LCA 

Y 

332 


120 

01 


LCA 




121 

01 


LCA 




122 

01 


W 




123 

08 

SKREC 

MCW 

( U3 

REC 

R 

124 

08 


B 

SKPON 

REC 

P 

125 

05 

PNEX 

CU 

( U3 


B 

126 

04 


8 

PEXIT 



127 

08 

SKPON 

B 

PNEX 

REC 

+ 1 T 

128 

04 


a 

SKREC 




- COMMENTS 


SUBROUTINE *NXREC* 

CONSTRUCTS AND PRINTS LINE OF Y GRID-LABELS 
STORE RETURN AT *ENDNR*. 

CLEAR PRINT AREA. 

MOVE YSTART TO TEMPI. 

CLEAR INDEX A. 

CLEAR INDEX B. IF FIRST 5 POSITIONS OF YSTART (NOW IN TEMPI) ARE 
NOT ZERO, GO TO *EDYST*. IF THEY ARE ZERO, MOVE LAST DIGIT 
(WITHOUT SIGN), AND IF IT IS ALSO ZERO, GO TO *YZROK* . 


LOAD Y EDIT WORD INTO TEMP2 AND EDIT TEMPI. 

LOOP TO LOCATE FIRST NON-BLANK POSITION IN YSTART (NOW IN TEMP2 ) • 
EXIT TO *NONB* WITH INDEX B = NUMBER OF BLANKS. 

SET LENGTH OF LABEL FIELD = 8. MOVE SIGN INTO BTH POSITION. 


THE EFFECT OF THIS LOOP IS TO BUILD A CONSTANT IN INOEX A, TO BE 
USED IN CENTERING A Y GRID-LABEL, AND A CONSTANT IN*TEN* TO 
BE USED IN *YZRCK* WHICH ADVANCES INDEX A REAOY FOR CENTER- 
ING THE NEXT LABEL. 


MOVE Y GRID-LABEL TO PRINT AREA USING CONSTANT IN INDEX A. 

INCREASE INDEX A BY CONSTANT IN *TEN*. 

RESTORE *TEN* = 10. 

ADD DY TO THIS Y GRID-LABEL TO GET NEXT LABEL. 

IF PRINT LINE IS NOT COMPLETED, GO BACK TO *NXTLAB*. WHEN IT IS 
COMPLETED, PRINT THE LINE OF LABELS AND EXIT THROUGH RETURN 
STORED AT *ENDNR*. 


ERROR EXIT. 

WRITE ERROR MESSAGE AND HALT THIS PLOT. 

PRINT A LINE. 

CLEAR PRINT AREA. 

MOVE CURRENT* VALUE OF Y INTO PRINT AREA. 

MOVE CURRENT VALUE OF ( X-LASTX ) , WH ICH IS NEGATIVE, INTO PRINT AREA. 
MOVE MESSAGE NAMED XOFORD INTO PRINT AREA. 

PRINT ERROR MESSAGE. 

READ A TAPE RECORD. 

IF FIRST CHARACTER IS P, GO TO *SKPUN*. 

IF NOT, BACKSPACE TAPE ONE RECORD AND RETURN TO TTP THROUGH 
♦PEXIT*. 

IF SECOND CHARACTER IS T, GO TO *PNEX*. 

IF NOT, GO TO *SKREC*. 


1405 

1409 

1413 

1414 
1421 

1428 

1435 

1442 

1447 

1454 

1462 

1469 

1476 

1484 

1491 

1495 

1499 

1506 

1514 

1521 

1528 

1535 

1539 

1546 

1553 

1560 

1567 

1574 

1579 

1583 


1587 

1588 

1592 

1593 
1600 
1601 
1602 

1603 

1611 

1619 

1624 

1628 

1636 



PLOT 


W 


NU. CC 

LABEL OP 

A-ADD. 

B-ADD. 

D 

- 

- 

- 

- 

- 

- 

- COMMENTS - 

------ 

LOCATION 













PROCESS INFO RECORDS 
(ENTRY AT ^FORMAT*) 



129 

04 

MULCV 

SW 

MCBRA 





SET SWITCH MCBRA FOR MULTIPLE CURVES. 

1640 

130 

04 


B 

MCBACK 





GQ TO *MCBACK*. 

1644 

131 

07 

FORMAT 

MCM 

REC 

+ 2 

IREC 



MOVE INFO RECORD (OMITTING 1ST 2 CHARACTERS -PI-) TO IREC. 

1648 

132 

07 


CW 

MCBRA 


LSTLAB 



SET SWITCH MCBRA FOR SINGLE CURVES AND SWITCH LSTLAB FOR FIRST- 

1655 

133 

08 


B 

MULCV 


PLCHAR 


- 

LEGEND. IF PLOT CHARACTER IS (-), GO TO *MULCV* AND RETURN. 

1662 

1 34 

07 

MCBACK 

SBR 

A 


0 



CLEAR INDEX A. 

1670 

135 

07 


MCW 

PXSH 


A 



CONSTRUCT THE EDIT WORDS , EDX AND EDY, FOR X AND Y GRID-LABELS, 

1677 

136 

07 


MCW 

PERZ 


EDX -6 

A 


USING DECIMAL-POINT-SHIFT CODES IN INDEX A TO POSITION THE 

1684 

137 

07 


MCW 

PYSH 


A 



POINTS. 

1691 

138 

07 


MCW 

PERZ 


EDY -6 

A 



1698 

139 

07 


M2 

XSTART 

-8 

XSTART 



MOVE SIGNS OF XSTART AND YSTART FROM HIGH TO LOW ORDER POSITIONS. 

1705 

140 

07 


MZ 

YSTART 

-8 

YSTART 




1712 

141 

07 


MCW 

ZTST 

-4 

XSTART-6 



REPLACE THE DECIMAL POINT IN XSTART WITH A ZERO. 

1719 

142 

07 


MCW 

DY 


DYT -1 



MULTIPLY DY BY 10 BY MOVING INTO AND OUT OF DYT. 

1726 

143 

07 


MCW 

OYT 


DY 




1733 

144 

07 


MZ 

DX 

-8 

DX 



MOVE SIGNS OF DX AND DY FROM HIGH TO LOW ORDER POSITIONS. 

1740 

145 

07 


MZ 

DY 

-8 

DY 




1747 

146 

07 

DXPL 

SBR 

A 


0 



CLEAR INDEX A. 

1754 

147 

07 


A 

ZTST 

+ 3 

YGD 



REPLACE LEADING BLANKS IN XGD AND YGD WITH ZEROES BY CHAINED 

1761 

148 

01 


A 






ADDING OF A SIGNED CONSTANT ZERO. 

1768 

149 

07 


C 

XGD 


A 



TEST X GRID-LINE FREQUENCY. 

1769 

150 

05 


B 

GRDX 




/ 

IF NON-ZERO, GO TO *GRDX*. 

1776 

151 

07 


MCW 

N INS 


XGD 



IF ZERO, REPLACE BY 999. 

1781 

152 

07 

GROX 

C 

YGD 


A 



TEST Y GRID-LINE FREQUENCY. 

1788 

153 

05 


B 

GRDY 




/ 

IF NON-ZERO, GO TO *GRDY*. 

1795 

154 

07 


MCW 

NINS 


YGD 



IF ZERO, REPLACE BY 999. 

1800 

155 

07 

GRDY 

A 

FM1 


A 



ADD LEFT-SIDE-FIELD LENGTH TO INDEX A. 

1807 

156 

07 


MN 

A 


PLOTA +3 



MODIFY ADDRESS OF *PLOTA*. 

1814 

157 

01 


MN 







1821 

158 

07 


C 

A 


ZTST -4 



TEST LEFT-SIDE-FIELD LENGTH. 

1822 

159 

05 


B 

MVNOP 




S 

IF ZERO, GO TO *MVNOP*. 

1829 

160 

07 


MCW 

« 

-6 

PLOTA 



IF NON-ZERO, MAKE *PLOTA* A MOVE COMMAND. 

1834 

161 

07 

MNVR 

SBR 

A 


3 

A 


ADD 3 TO INDEX A. USE NEW TOTAL TO MODIFY *XADU*. 

1841 

162 

07 


MN 

A 


XADD +3 




1848 

163 

01 


MN 







1855 

164 

07 


SBR 

A 


3 

A 


AGAIN ADD 3 TO INDEX A AND MODIFY ADDRESS OF *YADD*. 

1856 

165 

07 


MN 

A 


YADD +3 




1863 

166 

01 


MN 







1870 

167 

07 


SBR 

TWONE 


0 

A 


STORE INDEX A AT TWONE (6 + LENGTH OF LEFT-SIDE-FIELD). 

1871 

168 

05 


CC 

NXREC 




J 

DOUBLE SPACE AND GO TO *NXREC* TO PRINT A LINE OF Y GRID-LABELS. 

1878 

169 

04 


SW 

PLZERO 





SET A WORD-MARK IN PRINT AREA AT FIRST Y GRID-LINE (PLZERO). CLEAR 

1883 

170 

01 


cs 






ALL POSITIONS TO THE LEFT AND PUT THE X GRID-LINE CHARACTER 

1887 

171 

07 


MCW 

MINUS 


335 



(MINUS SIGN) IN ALL POSITIONS TO THE RIGHT. 

1888 

172 

07 


SBR 

MODXGD+3 

TEN 



FORCE FIRST X GRID-LINE AND LABEL, AND Y GRID-CHARACTERS (EVERY 10 

1895 

173 

07 


MCW 

TEN 


XL8LTS 



POSITIONS), BY PRE-SETTING * MOO XGD* , XLBLTS, AND YGDTST. 

1902 

174 

07 


MCW 

XGD 


YGDTST 




1909 

175 

04 


B 

START 





RETURN TO TTP THROUGH *START*. 

1916 



PLOT 


CD 

O 


NO. CC 

LABEL OP 

A-ADD. 

B-ADD. 

D 

— 






- COMMENTS - 




LOCATION 













PROCESS 

LEGENO 

RECORDS 



176 

01 

LABEL 

NOP 




LSTLAB=FIRST- LEGEND SWITCH. 

1920 

177 

04 

LSTLAB 

B 

LAB2 



IF NOT-FIRST LEGEND, GO TO *LAB2* • 

1921 








IF FIRST-LEGEND, ASK IF LINE WAITING IN PRINT AREA IS A 


178 

08 


BBE 

LABP 

331 

$ 

GRID-LINE. IF IT IS, GO TO *LABP* • 

1925 

179 

01 


W 




IF NOT, PRINT IT. 

1933 

180 

04 


CS 

332 



CLEAR PRINT AREA. 

1934 

181 

01 


cs 





1938 

182 

07 

CPYGD 

A 

DX 

XSTART 


INCREASE XSTART (VALUE USED FOR X GRID-LABELS) BY DX. 

1939 

183 

07 


A 

ONE 

XLBLTS 


STEP X GRID-LABEL COUNTER (XLBLTS) BY 1. 

1946 

184 

07 


C 

XLBLTS 

TEN 


TEST X GRID-LABEL COUNTER TO SEE IF IT IS TIME FOR 

1953 








A LABEL (EVERY TENTH LINE). 


185 

05 


B 

ENDGD 


S 

IF TIME TO PRINT AN X GRID-LABEL, GO TO *ENDGD*. 

i960 

186 

04 


B 

VTGRID 



IF NOT, GO TO SUBROUTINE * VTGRID* FOR Y GRID-LINES. 

1965 

187 

04 


VI 

CPYGD 



CONTINUE PRINTING GRID-LINES UNTIL IT IS TIME FOR NEXT X 

1969 








GRID-LABEL. 


188 

04 

ENDGD 

sw 

PL2ER0 



SET A WORD MARK IN PRINT AREA AT FIRST Y GRID-LINE (PLZERO). PUT 

1973 

189 

07 


MCW 

MINUS 

335 


X GRID-LINE CHARACTER (-> IN ALL POSITIONS TO THE RIGHT. 

1977 

190 

07 


SBR 

MODXGD+3 

TEN 


SET A-ADD. OF *M0DXGD*=10. 

1984 

191 

04 


B 

VTGRID 



GO TO SUBROUTINE *VTGRID* TO PUT Y GRID-LINE CHARACTERS INTO 

1991 








EVERY TENTH POSITION IN LAST X GRID-LINE. 


192 

04 

LABP 

B 

XGL 



GO TO SUBROUTINE *XGL* FOR X GR ID-LABEL. PR INT LAST X GRID-LINE AND 

1995 

193 

04 


W 

NXREC 



GO TO SUBROUTINE *NXREC* FOR LINE OF Y .GRID-LABELS. DOUBLE 

1999 








SPACE. 


194 

02 


CC 



L 


2003 

195 

04 

LAB2 

cs 

332 



CLEAR PRINT AREA. 

2005 

196 

01 


cs 





2009 

197 

07 


MCM 

REC +2 

201 


MOVE LEGEND RECORD TO PRINT AREA. PRINT AND RETURN TO TTP VIA 

2010 

198 

04 


w 

PEXIT 


• 

•PEXIT*. 

2017 








PATCH FOR PROCESS INFO 


199 

200 

07 MVNOP 
04 

MCW 

B 

FSTREC-1 

MNVR 

PLOTA 

MAKE *PLOTA* A NOP COMMANO. 
GO TO *MNVR* • 


2021 

2028 


PLOT 


NO. CC 

LABEL OP 

A— ADO. 

B-ADD. 

0 

— 

- 

— 

- COMMENTS - 

— _ 

LOCATION 










SUBROUTINE »VT GRID* 










INSERTS Y GRID-L INE CHARACTERS 

ID INTO PRINT AREA. 



201 

04 

VTGRID 

SBR 

YGXIT +3 



STORE RETURN AT *YGXIT*. 

2032 

202 

07 


SBR 

B 

0 


CLEAR INDEX B 

2036 

203 

07 

XGRIO 

MCW 

VERT 

PLZERO 

B 

MOVE Y GRID-LINE CHARACTER INTO PRINT AREA USING INDEX B. ADD 

2043 

204 

07 

MOOXGD 

A 

XGD 

B 


TO INDEX B THE QUANTITY SPECIFIED IN A-ADD. OF *MOOXGU*. 

2050 

205 

07 


C 

ONETEN 

B 


CONTINUE UNTIL B = 105. 

2057 

206 

05 

NOXGDA 

B 

XGRID 


T 


2064 

207 

07 


SBR 

MODXGD+3 

YGD 


RESET A-ADD. OF #MODXGD* TO Y GRID-LINE FREQUENCY. 

2069 

208 

04 

YGXIT 

B 

• 



RETURN 

2076 









SUBROUTINE *XGL* 









CONSTRUCTS AND EOITS X GRID-LABEL USING XSTART 


209 

04 

XGL 

SBR 

RETXL +3 



STORE RETURN 

2080 

210 

07 


C 

XSTART-1 

ZTST -2 


IF FIRST SIX DIGITS OF X START ARE NOT ZERO, GO TO *XGLA* • 

2084 

211 

05 


B 

XGLA 


/ 


2091 

212 

07 

XZRTST 

MN 

XSTART 

PLZERO— 3 


MOVE 7TH DIGIT (WITHOUT SIGN). IF IT IS ZERO, GO T0*RETXL*-7. 

2096 

213 

08 


B 

RETXL -7 

PLZERO-3 

0 


2103 

214 

07 

XGLA 

LCA 

EDX 

PLZERO-2 


LOAD X EDIT WORD INTO PRINT AREA AND EDIT XSTART. 

2111 

215 

07 


MCE 

XSTART 

PLZERO-2 



2118 

216 

07 


SBR 

XLBLTS 

0 


SET X GRID-LABEL COUNTER = 0 AND RETURN. 

2125 

217 

04 

RETXL 

B 

• 



RETURN 

2132 


O') 


i 


PLOT 


NO. 

CC 

LABEL 

OP 

A-ADD. 


------ DATA DEFINITIONS ------ 

LOCATION 

218 

01 

PPGUP 

DC 

• 


SKIP- TO— NEW- PAGE SWITCH. A WORD MARK IS SET AT PPGUP BY PLOT EACH 
TIME IT IS CALLED, FOR LATER UoE BY T TP TO START NON-PLOT 
OUTPUT ON NEW PAGE. 

2136 

219 

02 

PERZ 

DCW 

• 

0. 

CONSTANT USED IN CONSTRUCTING THE tDlT WORDS EDX AND EDY. 

EDY, EDX, EDIT 

THE NEXT THREE FIELDS MUST BE ADJACENT AND IN 
THE ORDER SHOWN. (USED IN A CHAINED MOVE.) 

2138 

220 

09 

EDY 

DCW 

• 

- 

EDIT WORD FOR THE Y GRID-LABELS. 

2147 

221 

09 

EDX 

DCW 

* 

- 

EDIT WORD FOR THE X GRID-LABELS. 

2156 

222 

09 

EDIT 

DCW 

• 


CONSTANT USED TO INITIALIZE THE EOIT WORDS EDX AND EDY. 

XOFORD , X, Y 

THE NEXT THREE FIELDS MUST BE ADJACENT AND 
IN THE ORDER SHOWN. (USED IN CHAINED MOVE.) 

2165 

223 

15 

XOFQRD 

DCW 

• 

X OUT OF ORDER 

ERROR MESSAGE. 

2180 

224 

03 

X 

DCW 

• 

000 

INPUT - THE LINE COUNT FOR A POINT OR SET OF POINTS. 

2183 

225 

03 

Y 

DCW 

• 

000 

INPUT - THE PRINT POSITION COUNT FOR A POINT. 

2186 

226 

03 

NINS 

DCW 

* 

999 

CONSTANT REPLACES XGD OR YGD IF EITHER IS ZERO. 

YGDTST, LASTX, XLBLTS 

THE NEXT THREE FIELDS MUST BE ADJACENT AND 
IN THE ORDER SHOWN. (USED IN A CHAINED ADD) 

2189 

227 

03 

YGDTST 

DCW 

• 

+ 00 + 

THE X GRID-LINE COUNTER. IT IS STEPPED BY ONE WHENEVER A LINE IS 
PRINTED AND COMPARED TO XGD. WHENEVER THEY ARE EQUAL, AN X 
GRID-LINE IS PRINTED AND THE COUNTER IS RESET TO ZERO. 

2192 

228 

03 

LASTX 

DCW 

* 

+ 00 + 

THE LINE COUNTER. IT IS INITIALIZED TO ZERO WHEN THE FIRST DATA 
RECORD IS PROCESSED AND STEPPED BY ONE EACH TIME A LINE IS 
PRINTED. 

2195 

229 

03 

XLBLTS 

DCW 

* 

+ 00 + 

THE X GRID-LABEL COUNTER. IT IS SET TO TEN TO FORCE THE FIRST X 

LABEL TO PRINT, THEN RESET TO ZERO. EACH TIME A LINE IS PRINT- 
ED, IT IS INCREASED BY ONE AND COMPARED TO TEN TO DETERMINE 
IF A LABEL SHOULD BE PRINTED. WHEN EQUAL TO TEN A LABEL IS 
PRINTED AND THE COUNTER RESET TO ZERO. 

2198 

230 

07 

ZTST 

OCW 

• 

0000000 

CONSTANT ZERO. 

2205 

231 

09 

DYT 

DCW 

* 

000000000 

FIELD USED TO MULTIPLY DY BY 10 BY SHIFTING IT LEFT ONE POSITION. 

2214 

232 

03 

TEN 

DCW 

* 

010 

CONSTANT. ALSO USED AS A COUNTER IN CENTERING Y GRID-LABELS. 
VERT, , ONE 

THE NEXT THREE FIELDS MUST BE ADJACENT AND 
IN THE ORDER SHOWN. (USED IN A CHAINED ADD) 

2217 

233 

234 

01 

03 

VERT 

OCW 

OCW 

• 

* 

1 

001 

THE X GRID-LINE CHARACTER. 

2218 

2221 

235 

03 

ONE 

DCW 

* 

001 

CONSTANT. 

2224 

236 

03 

ONETEN 

OCW 

* 

105 

CONSTANT USED TO STOP PRINTING Y GRID-LINES. 

2227 

2 37 

01 

RM 

OCW 

* 

4= 

CONSTANT. IF Y IS GREATER THAN 102 THIS CHARACTER IS PRINTED IN 
PRINT POSITION 132. 

2228 

238 

07 

TEMPI 

OCW 

• 


TEMPORARY STORAGE USED IN COMPUTING Y GRID-LABELS. 

2235 

2 39 

03 

199 

OCW 

• 

19 I 

CONSTANT USED IN CENTERING THE Y GR I D-LABELS . ( COMPL EMENT OF I) 

2238 

240 

03 

YL AST 

DCW 

* 

103 

CONSTANT. IF Y IS GREATER THAN YL AS I A RECORD MARK (t ) IS PRINTED 
IN PRINT POSITION 132. 

2241 

241 

03 

TWONE 

OCW 

• 

000 

COMPUTED DURING THE PROCESSING OF THE INFO RECORD, THIS IS 6 + FM1 • 
IT IS USED TO LOCATE THE NEXT OATA BLOCK IN THE DATA RECORD 

2244 


FOR SINGLE CURVES AND THE FIRST PLOTTING CHARACTER IN THE 
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APPENDIX A 


SAMPLE PLOTS 

The plots in figures 1 to 12 on the folio-wing pages were produced by 
specially written FORTRAN programs to illustrate some of the possibilities of 
the OPP system and to provide sample tape records for the 1401 systems pro- 
grammer. The same 7090 output tape read by the 1401 when the plots were gene- 
rated was also read by a specially written tape-dump program to present the 
associated tape records exactly as they appear on the tape,, except for the in- 
sertion of blank lines between records for easier reading. 
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DATA USED FOR THIS PLOT 


X - COS Y AND Y = 01.06)1.5 



PT FIGURE I- • 

PT PLQTXY KODE = 0 

P I *33 Q.OOIOOO 0.000000-0.000020 0.000020 10 10 0 • 

PD 0 0 0 3 0 6 1 9 1 12 2 15 3 18 4 21 6 24 7 27 9 30 II 33 12 36 14 39 17 42 19 45 21 48 24 51 26 54 29 57 32 60 
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PL DATA USED FOR THIS PLOT X = COS Y ANO Y = 0(.06)1.5 


Figure 1. - This plot illustrates using PLOTXY with KODE = 0. The subroutine provides 10 x 10 grid-line spacing, chooses 
both scales, and uses an asterisk as the plotting character. The number of points (26.) was placed in P(l). 
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Figure 2. - This plot illustrates using option 1 - choice of plotting character. The desired plotting character. 
X, was placed in P(2). Option 4 - choice of y grid-line frequency - was also used. The desired frequency (5.) 
was placed in P(4). Notice that the frequency of y grid-labels is not affected. 
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Figure 3. - This plot illustrates using option 32 - choice of y scale. The scale chosen was specified to PLOTXY by 
placing KSY = 3. in P(9), FY = -500. in P(lO), and DY = 25. in P(ll). The maximum possible value of KSY was used. 
(Fig. 4 displays the effect of using a smaller value of KSY.) Also illustrated is the use of option 8 - choice of 
printing left-s ide fi eld. The message was plac ed in the P array; the characters Ibbbbbbl in P(12J through P(l5), 
the characters 10PTI0NI in P(l6); the characters fbFhlSbi in P(l7); etc. 
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Figure 4. - This plot illustrates using option 64 - choice of printing coordinates at left. Because of the limited 
field size allotted to this print-out, PliOTXY transforms the coordinates to preserve as many significant figures 
as possible, and indicates in a special heading what has been done. This plot also illustrates using option 32 - 
choice of y scale. The same starting- value and scale-factor were wanted as in figure 3, but the values KSY « 1. 
were placed in P(9), FY = -50000. in F(lO), and DY = 2500. in P(ll). Notice the corresponding change in the 
appearance of the y grid-labels. 
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Figure 5. - This plot illustrates using option 2 - choice of x grid-line frequency. The desired frequency (20.) 
was placed in P(3). Notice that the frequency of the x grid-labels is not affected. This plot also demonstrates 
that after the last point has been plotted, the grid is only continued to the next grid-label, not to the next 
grid-line. The result of using options 8 and 64 simultaneously should be compared with figures 3 and 4. 
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Figure 6. - This plot illustrates using PLOTMY - Variation III - NO DUP with KKK(l) = KODE =0. As in PLOTXY, the 
grid-line spacing 10 x 10 is supplied by the routine, as are the scales and the plotting characters. The number 
of curves (3) was placed in KKK(2), the number of points for the first curve in KKK(3) , for the second in KKK(5), 
and for the third in KKK(7 ) . Variation III was specified by placing 5. in P(l). Notice that several points fall 
on the grid-lines and compare with figure 7. 
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Figure 7. - This plot illustrates using options 2 and 4 to remove all but the first, last, and left-most grid-lines, 
by supplying the frequencies = 0 in P(3) and P(4). This effect is possible in PL0TXY also. This plot also illu- 
strates the use of option 64 with PLOTMY . The x ordinate only is printed, not both, as in PLOTXY (see fig. 4). 
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Figure 8. - This plot illustrates using PLOTMY - Variation I - DUP X with KKK(l) - KODE =0. A single set of x values 
was placed in the x array, and five sets of y values one after another in the y array. The number of curves was 
placed in KKK{2) and the number of x values in KKK(3). Variation I was specified by placing 1. in P{l). 
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FIGURE 9. PLOTHY DUPY KODE 
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Figure 9. - This plot illustrates using PLOTMY - Variation II - DUP Y with KKK(l) * KODE = 16. The x and y arrays used 
for figure 8 were interchanged and the x scale specified hy placing KSX » 2. In P(6) and DX » 25. in P(8). Otherwise, 
PLOTMY’ s choice of scale (which is based on the length of the x array, not on the number of points In the y array 
supplied in KKK(3) ) would have spread the plot over four pages. Variation II was specified by placing 3. in P(l). 
Compare with figure 8. 
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PL (SAME OATA AS IN FIGURE 8.) 


Figure 10. - This plot illustrates the use of PLOTXY for multiple curves. The data used for figure 8 were supplemented 
by repeating the set of x values in the x array, then SORTXY was called. Option 16 was used and the x scale specified 
by placing KSX =5. in P(6) and DX - 1. in P(8) to avoid spreading the plot over four pages. (PLOTXY would have 
chosen a scale based on the number of points (115.) supplied in P(l) . ) Compare with figure 8. 
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PL 


Figure 11. - When either PLOTXY or PLOTMY chooses a scale, the value zero does not necessarily fall on a grid-line. If 
it is desirable to see intersections at a glance, (in this case, with the x-axis), this is readily done. For this plot, 
PLOTMY - Variation III - NO DUP was used, and a faked curve generated whose x values were . 1(. 1)4.9, {the range of x 
was known) and whose y values were zero. Option 1 was used in order to choose a zero for the plotting character for 
the second curve; option 16 was used to confine the curve to one page ; option 4 was used with P{4) = 20. to lighten the 
appearance of the grid. 
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PL * 


Figure 12. - This plot simulates a plot made on semilog paper. The values in the y array were those that would be plotted 
on the linear scale, while those in the x array were the logarithms of the values that would be plotted on the logarithmic 
scale. The original values are displayed using option 8. 



APPENDIX B 


NOTES 


I. Machine Configuration 

A. The FORTRAN subroutines have no hardware requirements other than those 
required for FORTRAN II (see IV. )j however, they do rely on the Lewis 
Monitor Error Package (see III.-B.), which makes use of divide-by-zero - 
trap hardware. 

B. The SPS subroutine PLOT was written for a 4K-1401, modified 1 to alter 

address -validity and instruction-length checking, and with the following 
"Special Features" (ref. 2): ^ 

Print Storage 

Additional Print Control 

Sign Control Left 

Decimal Control 

Indexing 

Store Address Register 
Move Record 
High-Low-Equal Compare 
Column Binary 

II. Storage Requirements 

PLOTXY - 875 words 
PLOTMY - 1125 words 
PISTUG - 375 words 
SORTXY - 97 words 
SCALE - 118 words 
PLOT - 1924 characters 


1 IBM-RPQ F -93971. 
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Ill* Operating and Programming Conventions 


A. The Lewis Monitor System requires all output to be written on tape 6. 
WRITE OUTPUT TAPE 6 statement numbers are: 

PLOTXY - 238, 251, 410, 420, 430, 440 

PLOTMY - 306, 318, 432, 462 

PISTUG - 900 

B. The Lewis Monitor System includes the Lewis Monitor Error Package, 
which is responsive to overflows, major and minor underflows, division 
by zero, and the misuse of library subroutines, such as, square root of 
negatives, logarithm of zero (see XII., Error Write -Outs ) . 

C. The Lewis 1401 tape-to-printer program expects input records on 
tape 3. In subroutine PLOT, instruction *FXREC* reads tape 3, and 
*FNEX* backspaces tape 3. 

D* The Lewis 1401 tape-to-printer program loads the word marks for the 
three index registers; hence subroutine PLOT does not. 

E. The Lewis 1401 tape-to-printer program tests a switch, *\PPGUP*, set by 
subroutine PLOT at *EUTRY*+7, to control skipping to a new page. 

F. Subroutine PLOT moves records with a Move Record command; hence all 
TITLE, INFO, and LEGEM) records must be terminated with a record mark 
or a group mark word mark. 

G. If the printing of a plot is halted by machine error or operator in- 
tervention, the subroutine PLOT must be reloaded. Reloading will supply 
the necessary initialization that occurs when a plot ends normally. 

IV*. FORTRAN II 

All FORTRAN subroutines were written in 32K 709/7090 FORTRAN (ref. 3) . 

A. There are Boolean statements in: 


PLOTKY 

- 140, 

150, 

160, 

170, 

180, 

190, 

200 

PLomr 

- 103, 

104, 

105, 

106, 

120, 

150, 

154, 160, 164, 170, 180, 


230, 

240, 

322 






B. The standard FORTRAN library subroutines for LOG and EXP are used by 
PLOTXY, PL0TMY, and PISTUG. 



C. The standard FORTRAN library subroutine IOH, modified to write 132 
character records, is used by ELOTXY and PLOTMY. 
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D. The standard FORTRAN library subroutine IOH will not generate a record 
less than three words (18 BCD characters) long, and it will supply blanks 
if necessary to attain that length. When the last PD record written con- 
tains only one pointy and the left-side field length is zero ; the end-of- 
record test made by PLOT fails, and the blanks are treated as the line 
count (zero) for another point. Although all points will have been 
plotted correctly, the plot will terminate with an X OUT OF ORDER error 
message and will not be properly ended* This situation can occur in 
PLOTXY only when conditions (l), (4), and (5) are simultaneously true, 
and in PLOTMY only when (2), ( 3 ), (4), and (5) are simultaneously true: 

(1) The number of points to be plotted is one more than an integral 

multiple of 21. 

(2) The number of points to be plotted -is one more than an integral 

multiple of 16. 

(3) Either Variation II (DUPY) or III (NO DUP) is selected. 

(4) Option 8 is not being used. 

(5) Option 64 is not being used. 

These combinations of conditions are relatively rare; each has occurred 
only once in over a year of use. 

E. The simultaneous use of a symbol as a fixed-point variable and as the 
index of a loop within the list of an output statement is accepted and 
compiled correctly by both 704 FORTRAN and 32K 709/7090 FORTRAN. The 
simplest change for other versions of FORTRAN is the substitution of JJ 
for each appearance of J in the output statements 410, 420, 430, and 440, 
in PL0TXY . 

V. Perforation Test (PISTUG) 

When an x array is being scaled by PISTUG, the maximum number of lines 
permitted for the plot is computed by dividing the total number of points by 
35 (using integer arithmetic), adding 1, and multiplying the result by 55. 
(The constants 35 and 55 were determined empirically.) After the scaling 
parameter for the scale-factor has been computed, a check is made to deter- 
mine whether, on a two-page plot, the final x grid-line will fall within 12 
lines of the top of the second page. If it will, an arbitrary 35 is added 
to the maximum number of lines permitted, and the parameter is recomputed. 

VI. Range Test (PISTUG) 

The absolute value of the range (maximum element minus minimum element) 
of any array being scaled by PISTUG must be greater than 10““^, and the 
larger of the absolute values of these two elements must fall between 10~ 3 
and 10?. This restriction is dictated by the length of the grid-label field 
on the 1401 and is enforced by testing the size of KCDX in PISTUG (see 
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XXI* -D.). If scaling is done by the programmer, no check is made to 
enforce this requirement, and the programmer is responsible for avoiding 
a choice of scale that will lead to impossible grid-labels. 

VII. Range Assumptions (PLOT) 

A. PLOT assumes the line count, KX, is limited by 0 < KX < 999, and if 
PISTUG does the scaling, it is. If the programmer uses option 16 and 
chooses scaling parameters that lead to negative line counts, program 
may be overwritten, with unpredictable results. Values of KX greater 
than 999 also cause trouble, but any programmer with 15 pages of plot 
has trouble already. 

B. PLOT assumes the print position count, KX, is limited by 0 < KY < 102, 
and if PISTUG is doing the scaling, it is. PLOT tests KY, and any value 
of KY such that 103 ^ KY ^ 999 will cause a record mark to be printed in 
position 103. If KY is negative and a single digit, the point will be 
printed as if it had been positive and the programmer will never know it 
is wrong; larger negative values may overwrite program inviting cata- 
strophe. 

VIII. Range Assumptions (PISTUG) 

Theoretically PISTUG can scale an array of any length correctly, but if 
scaling an X array of more than 595 points, scaling parameters leading to 
line-counts greater than 999 may be chosen. This limitation is on the total 
number of points, and would not apply, for instance, to using PLOTMI - 
Variation I - DUP X for six curves of 200 points each. 

IX. Direction Assumptions (PLOTMT) 

PL0TMY assumes that the scale-factor in the x-direction is positive, 
and if PISTUG is doing the scaling, it will be. If the programmer uses 
option 16 and supplies a negative scaling parameter, LX, the plot will be 
incorrect. If the programmer really wants a decreasing x scale he should 
multiply the X array by -1. and manually insert negative signs and alter 
n -" to "+ TT on the x grid-labels as required. 

X. Duplicate Points 

Even though PL0TMY may supply two or more plotting characters to be 
printed in the same location, PLOT will print only the last character re- 
ceived. No over-printing is done. Similarly, of multiple left-side fields, 
only the last received by PLOT will be printed. 

XI. Multi -record Titles and Legends 

The most common programming error to date is the omission of the re- 
quired 2HPT or 2HPL for the second TITLE or LEGEND record from the FORMAT 
statements. Such omission also occurs when slashes (///) are used in the 


107 



FORMAT statement to get line spaces following, or between, printed lines. 
This does not work. Use 2HPT/2HPT/2HPT, etc. 

XII. Error Write-Outs 

A. PLOTXY and PLOTMY have error returns. If the value of either of the 
scaling parameters DX or DY is returned from PISTUG as zero, the routines 
branch immediately to RETURN and no IHFO or DATA records are written. 
However, this situation has yet to occur without some message having 
already been written during the execution of PISTUG. 

B. PLOTMY will write an error message 

NPTS IS MISSING OR WRONG IN THE KARRAY. IMIN = nnn 

if the search loop (statement number 424 to 430), which finds the correct 
plotting character, is unsuccessful. Information missing or misplaced 
in the K array, or "dimensioning" the K array less than 14, has caused 
this write-out to occur. 

C. The programming error of specifying DX (or DY) = 0 is caught by the 
Lewis Monitor Error Package, either as division by zero in statement 
number 324 (or 334) in PLOTXY, and in 394 (or 410) in PLOTMY j or, if 

FX (or FY) is zero (this is permissible) also, DX (or DY) = 0 is caught 
earlier as LOGF(O.O) at 187+4 (or 199+4) in PLOTXY, or at 238 in PLOTMY. 

D. PISTUG writes an error message if an array fails the grid-label-size 
test made on KCDX, 

PLOT HALTED BY ARRAY NUMBER n VALUES OUT OF RANGE. KCD=mmmm 

and branches to RETURN with DX (or DY) = 0, which prevents further calcu- 
lation in PLOTXY or PLOTMX. 

E. The programming error of supplying an array (to be scaled by PISTUG) 
whose elements are all equal (usually zero, from using a wrong array 
name) is caught by the Lewis Monitor Error Package as LOGF(O.O) at state- 
ment number 122 in PISTUG. Sending valid data and forgetting to supply 
the number of points has exactly the same effect, as does supplying the 
number of points in fixed point rather than in floating point. If the 
programmer chooses both his scales and neglects to supply the number of 
points, there will be no error write-outs, but his plot will consist en- 
tirely of one horizontal grid-line correctly labeled. 

F. PLOT writes the error message 

X OUT OF ORDER mmmnnn 

whenever the value of a difference it forms (by subtracting the current 
total in the line counter from the current value of X) is negative. This 
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will always occur if PLOTXY has been called with the elements of the 
X array not in sequence. Also written out are this difference (mmmm) 
and the corresponding value of Y(nnn) . PLOT then bypasses the remaining 
PD and PL records without using them, and control is returned to the 
tape -to-pr inter routine. 
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